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CN ' Abstract 

We study the two-point correlation function of a freely decaying scalar in Kraichnan's 

model of advection by a Gaussian random velocity field, stationary and white-noise in time 

but fractional Brownian in space with roughness exponent < C < 2, appropriate to the 

inertial-convective range of the scalar. We find all self-similar solutions, by transforming 

^~>' the scaling equation to Kummer's equation. It is shown that only those scaling solutions 
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with scalar energy decay exponent a < (d/7) -I- 1 are statistically realizable, where d is space 
dimension and 7 = 2 — C. An infinite sequence of invariants Jp, p = 0, 1,2, ... is pointed 
out, where Jq is Corrsin's integral invariant but the higher invariants appear to be new. 
Lj " We show that at least one of the invariants Jq or Ji must be nonzero (possibly infinite) 

3 . for realizable initial data. Initial data with a finite, nonzero invariant — the first being Jp — 

converge at long times to a scaling solution $p with a = ((i/7) + p, p = 0,1. The latter 
belong to an exceptional series of self-similar solutions with stretched-exponential decay in 
space. However, the domain of attraction includes many initial data with power-law decay. 
When the initial data has all invariants zero or infinite and also it exhibits power-law decay, 
then the solution converges at long times to a non-exceptional scaling solution with the same 
power-law decay. These results support a picture of a "two-scale" decay with breakdown of 
self-similarity for a range of exponents {d + 7)/7 < a < (d + 2)/7, analogous to what has 
recently been found in decay of Burgers turbulence. 
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1 Introduction 

A key hypothesis in the theory of turbulence decay is that the energy spectrum has asymptot- 
ically a self-similar form 

E{k,t)=v'^{t)l{t)F{ke{t)) (1.1) 

where i{t) is a suitable length-scale and v{t) a velocity scale. Equivalently, the hypothesis may 
be made for the longitudinal velocity two-point correlation BLL{r,t) := fifj{vi{x)vj{x + r)), 
that it obey: 

BLLir,t)=v\t)fir/i{t)). (1.2) 



Historically, hypotheses ( |1.1| ) and ( |1.2[ ) were first proposed for freely decaying turbulence in 
1938 by von Karman and Howarth |l|. They referred to the Ansdtze (1.1) and (1.2) as "self- 



preservation hypotheses," since the shape of the spectrum and correlation function are thereby 
preserved in the decay process. Detailed discussion of such self-preservation hypotheses is to 
be found in |^, Chapter 16. Since the early work in turbulence, corresponding hypotheses 
have been proposed for many other nonequilibrium processes, e.g. surface growth [^ and 
phase-ordering dynamics ^. In those fields the assumption is usually called dynamic self- 
sim^ilarity or dynam,ic scaling. Essentially, the hypothesis amounts to the statement that there 
is only one relevant length-scale in the decay process. For example, in turbulence decay this is 
plausibly the integral length-scale L[t) := -g — 5__^ J^ dj- BLL{r, t). Althoughly widely employed 
in nonequilibrium physics, the validity of such self-similarity hypotheses is still debated and their 
foundations poorly understood. This is particularly true when the random initial data of the 
system exhibit long-range power-law correlations. 

Recently, the validity of the self-similarity has been examined in a soluble model, the de- 
caying Burgers turbulence ||5[. Those authors solved exactly for the two-point correlations and 
energy spectra of the one-dimensional Burgers equation with initial energy spectra exhibiting 
a low-wavenumber power-law form, E{k,tQ) ~ Ak"^ for kLo <^ 1, and thus a power-law decay 
in the spatial velocity correlation function, B{r,tQ) ~ A'r^^"'^^' for r ^ Lq (when n is not 



a positive, even integer). What was discovered by the authors of [^] was that the hypothesis 
of dynamic self-similarity was violated when 1 < n < 2. Instead, a new length-scale L*(t) 
developed dynamically which was much larger asymptotically than the integral length-scale 
L(t). The new length-scale was characterized by the property that the initial low-wavenumber 
power-law spectrum was preserved only for kL^:{t) ^ 1. This preservation was traditionally 
believed to hold for all kL{t) ^ 1, which has been called the principle of permanence of large- 
eddies [^. The development of two distinct length-scales had important consequences for the 
decay process. For example, the rate of decay was found to be different than that predicted by 
traditional phenomenological theory. 

Another exactly soluble turbulence model is available, a model of a turbulently advected 
scalar proposed by R. H. Kraichnan in 1968 |Q]. This model corresponds to a stochastic 
advection-diffusion equation 

dt9{r, t) + (v(r, t)- V)0(r, t) = k A 0(r, t), (1.3) 

in which the advecting field v(r,t) is a "synthetic turbulence", specifically, a Gaussian random 
velocity field with zero mean and covariance 

{vi{T,t)vj{r',t')) = A,(r - r')5{t - t') (1.4) 

which is white-noise in time. The remarkable feature of the model which Kraichnan discovered 
is that there is no closure problem. In particular, Kraichnan showed that the 2-point correlation 
function 0(r, t) := {9{r,t)9{0,t)) obeys the following closed equation in homogeneous scalar 
decay: 

dtG{r,t) = [A,(0) - Ai(r)]Vr,Vr,e(r,t) + 2k A e(r,t). (1.5) 

Recently the study of the Kraichnan model has undergone a renaissance, impelled by the 
observation [^ that the A^th-order statistical correlations for N > 2 should exhibit "anomalous 
scaling", not predicted by naive dimensional analysis as in Kolmogorov's 1941 theory. In 
particular, attention has been focused on the inertial-convective range of the model, in which 



the molecular diffusivity k — > and the velocity covariance in space has a power-law form 
Dij{r)^D,,{0)-Di-r'^ 



1 
d-1 



^ii + 7 I ^ij 



OijY, (1.6) 



for < C < 2 and r <^ L, where the latter length-scale is the velocity integral scale. In 
this relation statistical isotropy as well as homogeneity has been assumed. The formula (1.6) 



mimicks the situation in a real turbulent scalar decay when the scalar is spectrally supported 
on length-scales Lg much smaller than the velocity integral scale L but yet much larger than the 
dissipation length-scale rjd set by the molecular diffusivity k. In this situation, the calculation 
of the anomalous exponents has proved to be possible analytically by perturbation expansion 
in three regimes: small space Holder exponent H = ^ corresponding to a "rough" velocity field 
1^, H^, small exponent 1 — H coresponding to expansion about a "smooth" velocity field or 
so-called Batchelor limit [0, |T^ , and expansion in ^ with d the space dimension [^, |l^ . 

Our interest here is instead to study the Kraichnan model as a soluble test case of turbulent 
decay. Passive scalars undergo a turbulent decay which is similar in many respects to that 
of the velocity itself. This is very well described in [p!5| ]. The decrease in the scalar energy 
or intensity E0{t) := ^(^^(0, t)) is by a process of progressive degradation of the scalar at 
the higher wavenumbers. In the course of this decay process, the scalar integral length-scale 
LQ{t) = Q7g-jT /q°° dr 0(r, t) grows as the spectral support of the scalar is shifted progressively 
to lower wavenumbers. On dimensional grounds, this growth of the length-scale is governed by 

1 dLg{t) ^ ^ ^_^ 



a DiL-^t), (1.7) 



Le{t) dt 

with 7 := 2 — C, which leads to the relation Lg{t) oc {Di{t — to))"^ • As has been empha- 
sized in [1^, this can be thought of as a "Richardson diffusion" of the scalar integral length- 
scale up through the velocity inertial range (which is here taken to be statistically stationary). 
This growth law for the scalar length-scale can be converted into an energy decay law under 
two additional assumptions. First, if one considers initial scalar spectra of a power-law form 
Eg{k,to) oc Ak'^ for low- wavenumbers kLg{to) <C 1, then the hypothesis of permanence of 



large eddies would imply that this low-wavenumber spectrum persists with a time-independent 
constant A for kL*g{t) ^ 1, if n < d + 1. Second, the hypothesis of dynamic self-similarity 
would imply that there is only one relevant length-scale, the scalar integral length Lg{t), so 
that Lg(t) = Lg{t) up to a constant factor. Under these two assumptions, the scalar energy 
may be estimated to order of magnitude by integrating over the low-wavenumber range up to 
Lg (t), with the result 

Ee{t) oc ^V(t) oc A{Diit - t^T^I^ (1.8) 

for p = n + 1. Thus, the decay rate is non-universal, and depends upon the low-wavenumber 
spectral exponent. 

The main aim of this work is to examine the decay problem in the Kraichnan model, 
to investigate the universality of the "two-scale" phenomenon discovered in j^. Because the 
standard phenomenology is common to both velocity and passive scalar decay, we may use the 
Kraichnan model as a source of insight. There seems to have been less work on the decay 
problem in the Kraichnan model than on the statistical steady-state and most of this in the 
Batchelor limit C = 2, the so-called viscous-convective range. In addition to the early work 
of Kraichnan [|^, the decay of the passive scalar in the Batchelor limit has been recently 
examined by Son [^]. Our work here will be devoted instead to the inertial-convective range of 
the scalar in which the velocity correlator is given by (|1.6| ) with < C < 2. A few preliminary 
investigations on this problem have been reported in a very recent review article of Majda and 
Kramer [^], Section 4.2, but no exhaustive study seems yet to have been made. Our interest is 
thus in the low-order N = 2 correlator, rather than in the higher-order statistics. In the same 
inertial-convective range considered above, with the assumptions of homogeneity and isotropy, 
equation ( |1.5D for the 2-point correlator simplifies to 



dMr,t) ""' ^ 
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We shall be particularly interested in the issue of self-similarity of the decay process. In fact, 
one of our main results will be an analytical construction and complete classification of all self- 
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similar decay solutions of equation (1.9), along with an analysis of their domains of attraction. 
The reader should note that in equation (^]^) the limit has been taken of vanishing molecular 
diffusivity. We have done so in order to focus on the turbulent dissipation of the scalar, which 
leads to a decaying scalar energy Eg{t) — > as t ^ oo even in the limit k — > 0. This corresponds 
to a famous conjecture on the three-dimensional energy cascade put forth by Onsager p!9| ], who 
proposed that the limiting turbulent ensemble in the limit of vanishing viscosity i/ — > should 
consist of realizations of the inviscid Euler equations which dissipate energy. Of course, these 
must be weak or distributional solutions, not classical solutions. Onsager coined the term "ideal 
turbulence" for this limiting dissipative ensemble governed by the ideal fluid equations. More 
recently, this ideal mechanism of dissipation has been called the "dissipative anomaly," since 
Polyakov pointed out a close similarity to conservation-law anomalies in quantum field theory 
This specifically turbulent mechanism of dissipation is well-illustrated by the Kraichnan 
model. Since the operator on the RHS of (|1.9| ) is homogeneous degree —7, it follows using 
Eg{t) = l@{0,t) that 



dEg Di d 



dt ^ ' 2r^-i dr 
when the 2nd-order structure function 



.^«-'f(M) 



= (1.10) 

r=0 



52(r,t) := ([^(r,t)-0(O,t)]2) 

= 2[e(0, t) - e(r, t)\ ~ Cr^, (1.11) 

with ^ > 7. Thus, some critical degree of singularity is required for turbulent dissipation. Our 
self-similar decay solutions — which are exact solutions of the zero-diffusion Kraichnan equations 
— explicitly illustrate this ideal dissipation mechanism. The implications for the theory of weak 
solutions of the hyperbolic stochastic PDE, the k — > limit of (|l.3| ), will be discussed in a 
forthcoming work [pi|]. 



2 Background Material 

(2. 1) Review of the Kraichnan Model 

The model we consider is the stochastic partial differential equation 

de{r, t) = KA e{r, t)dt - (W(r, dt)- V)0(r, t), (2.1) 

where W(-,t) is a Wiener process in the function space C(M'^,M''), with covariance function 

{Wi{v,t)Wj{r',t')) = Dij{r-r')tAt'. (2.2) 

The stochastic PDE is interpreted in the Stratonovich sense. See ^2|, Chapter 6, and ||2^ for 
a more detailed discussion of the mathematical foundations. The spatial covariance matrix D 
we consider is defined by the Fourier integral 

where < C < 2. The constant ko is an infrared cutofi' for the velocity field, proportional to 
the inverse velocity integral length /cq oc L^^. -Pij(k) is the projection in M'^ onto the subspace 
perpendicular to k. Thus ( |2.3| ) automatically defines a suitable positive-definite, symmetric 
matrix-valued function, divergence-free in each index. We have made the choice ( p.3| ) just for 
specificity. In fact, any velocity covariance with the properties discussed next would suffice. 
The matrix -Djj(r) can be written as -Djj(r) = Pij{'Vr)K{r), or as 

Aj(r) = K{r)6,j + d^djH{r), (2.4) 

where the function K{r) is defined by the integral 

Kir) = D [d'^k ""' ' ,,,, (2.5) 

and H{r) is given by the (for d = 2, principal part) integral 

1 e*^*"" 



H{r) = D I d^k^ ^ ,,^,w, , (2.6) 



so that — A H = K. The scalar function K{r) is essentiahy just the standard Bessel potential 
kernel |24|, and may thus be expressed in terms of a modified Bessel function: 



Kir) = D ^^^^ . {korfl^K^,2{kor). (2.7) 



2 

The Hessian matrix didjH{r) of the function H of magnitude r = |r| alone is 

didjH{r) = 6ijJ{r) + ?i?j ■ r^ir), (2.8) 

with J{r) = H'{r)/r and r = r/r. However, because Tr(V(X' V-ff) = —K, a Cauchy-Euler 

equation follows for J[r): 

r^(r) + d-J(r)=-K(r). (2.9) 

dr 

Due to the rapid decay of its Fourier transform, the function J(r) is continuous. Thus, the 

relevant solution is found to be 

J{r) = -r-'^ [ p'^-^K{p)dp. (2.10) 

Jo 

in terms of K{r). Using this expression for J(r), along with Eq.( |2.8| ), we thus find 

Dij{r) = {K{r) + J{r))5^j - {K{r) + d ■ J{r))?i?j, (2.11) 

which gives Dij as an explicit linear functional of K. Cf. [Q], equations (14.1), (14. 3). 

We are interested to consider the model in the range of length-scales r <C L. We require 
some asymptotic expressions in that range: 

K{r) =Ko- Kir'^ + O (fcgr^) , (2.12) 

with 

r (i\ 7r'^/2 



K, = D^±^.k-^ (2.13) 

and 

Ki = D ^-^^ ^. (2.14) 

2C.c-rte' 



Also, 



with 



and 



Dij{r) = Do6ij - Di ■ r< ■ 



^'^^d^ 



nrj 



+ 0{k, 



2„2^ 



Dn 



d-1 



-Kn 



d + C 



(2.15) 



(2.16) 



(2.17) 



The first equation is derived by means of the known Frobenius series expansion for the modified 
Bessel functions (e.g. |25], section 7.2.2, equations (12), (13)). These give 



z^K^{z) = l^-^^^^z^^ + 0(z-^ 



■)\-v 



■ 21+^ 



(2.18) 



From this expansion for K^{z) and from the representation ( p.7D for K{r) we obtain the asymp- 
totic expression (2.12). Next, using ( p. 11 ), we observe that, if K{r) has a power-law form, 
K{r) = Kr^, then it is easy to calculate that 



Dij{r) = Kr^ 



d-1 



Sii + 



e 



- {6ij - ?i?j) 



(2.19) 



Since Dij{r) is linearly related to K{r), we may apply this formula to the first two terms in the 
expansion ( 2.12| ), taking first ^ = and then ^ = C- This yields the second expansion ( |2.15| ). 

As already shown by Kraichnan 0, the equation for the 2-point correlation functions in the 
model (|2.lD is closed. This was subsequently generalized to the A^-point correlations |26, 27 1. 



By now, these equations have been derived by several arguments and in many places, e.g. |28|. 
Therefore, we shall give no derivation here. However, we make a few comments on the physical 
interpretation in the case N = 2. Since Dij{0) = D^dij, the two terms from Dij{0) — Dij{r) 
may be treated separately, with the result that ( |1.5D may be written as 

(2.20) 



dtQir,t) = (Z)o + 2k) AG(r,t) - Aj(r)Vr,Vr,e(r,t). 



Hence, we see that the first term gives essentially an augmentation of the molecular diffusivity, 
i.e. it produces an "eddy diffusivity" Keddy = ^Dq. In fact, it is the same eddy diffusivity which 
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appears in the A^ = 1 equation [^8|. The second term, as we discuss in more detail below, 
represents additional triadic interactions between one velocity mode and two scalar modes. 
From ( 2.13| )-( |2.17] ) we see that both of these terms are separately infrared divergent, in the 



limit ko -^ 0, but that this infrared divergence cancels in the equation ( |1.5| ) for 0(r, t). Since 
we are only interested in the inertial-convective range behavior of the scalar, when Lg{t) <^ L, 
it is convenient for us to take the limit /cq -^ 0. This has been done in deriving ( |1.9D , which 
follows easily from (1.5) and (|2.15| ). 



We shall also need below the equation for the spectral scalar energy transfer. We introduce 
the Fourier transform 



G(k,t):=^^y e(r,t)e-''-dr. (2.21) 

In terms of this, the scalar energy spectrum is 

Ee{k,t):=^iOd-ik''-^eik,t), (2.22) 

where u>d-i = 2tt ''^ /T (2) is the (d— l)-dimensional measure of the unit sphere in d-dimensions 
and, if the scalar statistics are not isotropic, Q{k,t) is a spherical average. It is straightforward 
to Fourier transform (|1.5|), with the result 



dte(k,t) = -kikj I (ft\D,j{ci) e(k,t)-e(k-q,t) - 2Kk^%{\i,t) 

= -{Do + 2K)k'^e{k, t)+ I A (k"^D(q)k) • G(k - q, t), (2.23) 

where Dij{q) = K{q)Pij{q). Notice that this last d x d matrix is positive semidefinite. Also, 
for each k, 0(k, t) > 0, as a statistical realizability requirement. Thus, we can now see that 
in the spectral representation for each wavevector k, the first L'-term is always negative and 
represents a "loss" term, while the second D-teiui is always positive and gives a "gain" term. 
The first is an "eddy diffusivity" effect, as we have already discussed. The second can be seen 
to result from triadic interactions of a velocity mode with wavevector q and two scalar modes 
with wavevectors k — q and k. It is easy to derive from this an expression for the transfer term 
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T0{k,t) in the spectral energy balance equation 

dtEg{k,t) =Tg{k,t) -2Kk'^Ee{k,t). (2.24) 

There are "loss" and "gain" terms which are both infrared divergent in the limits /cq ~^ 0, but 
the infrared divergence cancels exactly in the equation for Eg{k,t). 
(2.2) Phenomenology of Turbulent Scalar Decay 

As discussed briefly in the Introduction, the standard phenomenology of turbulent decay is 
built upon two fundamental hypotheses: the permanence of large eddies (PLE) and dynamic 
self- similarity (DSS). We shall review each of these topics here in turn, commenting upon the 
original motivations for these hypotheses and their dynamical justification (or not) within the 
Kraichnan model. 

First, we consider the permanence of large eddies. The motivation for this hypothesis 
lies in the phenomenon of spectral hacktransfer of energy. As was first shown by Proudman 
and Reid in a calculation with the quasinormal closure for decaying three-dimensional, ho- 



mogeneous turbulence |29|, the rate of change of the energy spectrum asymptotically at very 
low wavenumbers is dominated by a small but significant source of energy, which arises from 
nonlinear interactions of energy-range modes. Because the source-term is positive, and hence 
opposite in sign to the forward-cascading transfer through the inertial subrange, this phe- 
nomenon is called "backtransfer" . According to calculations within spectral closures — such as 
quasinormal closure or its more sophisticated descendants, such as eddy-damped quasinormal 
markovian (EDQNM) closure — the transfer rate T{k, i) is a power-law form B{t)k'^~^^ in d space- 
dimensions. See ]2^ and also Q, Sections 15.5-15.6. The same phenomenon for scalar transfer 
that Te(A;, t) ~ B0{t)k'^^^ at low k was subsequently found by Reid using again the quasinormal 
closure [Q. See also [|l5[ and lH, Section 19.4. These closure calculations lead one to expect 
that, if the initial scalar spectrum has a power-law form, Eg{k, tg) ~ Ak'^ for kLQ{tQ) <^ 1, then 



this state of affairs will be preserved in time for n < d + 1. Indeed, with the latter assumption, 

fto' 
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the initial spectrum will dominate the time-integral f^ ds Tg{k,s) ~ BQ{t)k'^^^ asymptotically 



for small enough k. Hence, one may expect that Eg{k,t) ~ Aok"" for kLg{t) <C 1 with a con- 
stant independent of time. This is the usual statement of the hypothesis of permanence of large 
eddies, in a spectral formulation. 

One may also formulate a permanence hypothesis in physical space. Thus, in the scalar 
case, one may suppose that, if 0(r, t) ~ A'gr"^ for r ^ L{t) at the initial time t = to, then this 
relation will persist with the same constant A'g at later times t > Iq, at least when p < d + 2. 
In the case of velocity correlations for Navier-Stokes turbulence, Proudman and Reid showed 
that pressure forces induce a long-range power-law r"*-"^"*"^' at any positive time, even if such 
correlations are not present initially. This is the exact physical-space analogue of the k'^'^^ 
spectral backtransfer. Thus, whenp < d + 2 the correlations present initially shall dominate for 
r sufficiently large. Similar arguments apply to the scalar correlations. Although the long-range 
pressure forces drop out of the expression for the first time-derivative of the scalar correlation, 
they appear in the second- and higher-order derivatives. Thus, a similar power-law decay is 
expected there. 

Although there is a formal correspondence between the spectral space and physical space 
formulations of the hypothesis of permanent large eddies, the two versions are not equivalent. 



as has been emphasized in [g]. Formally (e.g. see [31|, Ch.IX.6, Theorem 4) 



corresponds to 



with p = n + 1 and 



E0{k,t) ~ Aek"^ for kLe{t) < 1 (2.25) 



e(r,t)~^^r-P for r»L0(t), (2.26) 



A', = 2P ^y^_\^M e. (2.27) 



Even formally, one can see that there is a problem when p — d = 2m, an even, nonnegative 
integer, since then A'q = given by the above formula, when Ag is finite. In fact, this case 
corresponds to Q{k,t) ~ (/c^)™ for kLg(t) <^ 1, which is thus analytic at small k. Hence, one 



13 



may expect the decay of 0(r, t) in physical space to be faster than any power, consistent with 
the vanishing of A'g in ( 2.27| ). It is therefore quite possible that Eg{k,t) exhibits a power-law 



(with or without permanent coefficient) and that 0(r, t) has no power-law behavior whatsoever. 
A case in point is when the function 0(r, t) is integrable, with a nonzero integral, 

< K{t) := j^ J dr eir,t) < c^. (2.28) 

This is consistent with rapid decay in physical space faster than any power, e.g. exponential. 
However, the energy spectrum in this case exhibits a power-law at low wavenumber. In fact, 

Eg{k,t) ^ A{t)k'^-^ (2.29) 

with 

A{t):=^LOd-.i-K{t), (2.30) 

just using the definition of the energy spectrum. This is sometimes called an "equipartition 
spectrum", because it represents an average modal energy which is the same for each Fourier 
mode k. Not only is there a power-law spectrum, but, in fact, the coefficient A{t) is independent 
of time t, i.e. spectral PLE holds. This is true because K{t) is known to be a constant of the 
motion, called the Corrsin invariant. It was first derived for the true passive scalar by S. 
Corrsin in 1951 |32, ^], who employed the equation 



2 8 
dte{r,t) 



^'^-^ I B{r,t) + ^—{r,t) 
or 



(2.31) 



with 

Bir,t)■.= ?■{^rir,t)9ir,t)9iO,t)). (2.32) 

This equation plays the role of the von Karman-Howarth equation for the temperature field 
and it is analogous to the equation ( |1.9| ) in the Kraichnan model. It is not hard to see that 
this equation implies that dK{i)/dt = when B{r,t) = o [r~^ >Y See Q, Section 15.2. 
Furthermore, the quasinormal closure calculations suggest that the indirect effect of pressure 
forces in dtB{r, t) (i.e. in the higher-order time-derivatives) lead to a decay at least as B{r, t) = 

14 



O (r~' "'"^^), provided that the decay was no slower initially. See [^, Section 15.6. The induced 
power-law decay in the scalar correlation itself is Qir^i^ = O (r^'^^+^^V which is integrable and 
thus the Corrsin integral should still be finite. 

All of these arguments, which are derived under some closure assumptions for true passive 
scalars, are derived more directly and convincingly in the Kraichnan model. For example, if 
initially there is a power-law decay of correlations in space as in ( 2.26| ) in the Kraichnan model, 
then 

ate(r,to)~-I)ip(d + 7-p)^e^"^P+^^ for r>Le(to), (2.33) 



since the operator on the RHS of ( |l.g| ) is homogeneous degree —7. Since 7 > 0, this decay is 
faster, and it does not seem possible for the initial power-law to be upset at large enough r. 
Notice that in the Kraichnan model — unlike for true turbulent scalars — there is no action of 
pressure forces on the scalar correlators whatsoever. Hence, in the Kraichnan model, spatial 
PLE should hold for all p > and not just for < p < d + 2. In the spectral formulation, 
however, there is such a restriction. It follows from (p. 23) that in the limit as A; ^ 0, 



dtEeik, t) ~ -{Do + 2K)k'^Ee{k, t) + Be{t)k'^^^ (2-34) 

with 

Bg{t) := [ d'^q (k^D(q)k) • e(q, t) > (2.35) 

and k any unit vector. Thus, the "gain" term provides exactly the k'^^^ power-law backtransfer 
term expected at low wavenumber k. Hence, spectral PLE should hold when there is a power- 



law spectrum of the form ( 2.25 ) initially, with n < d + 1. Finally, the theorem on the Corrsin 
invariant is easily seen to follow by multiplying (|1.9| ) by r'^^^ and integrating over r. The result 
follows as long as 0(r, t) is o (r"^'^"'^') near r = and r = 00. The moral is that PLE, whether 
in physical or spectral space, should hold when the necessary power-laws are present initially. 
The second major hypothesis invoked in the traditional theory of turbulence decay is dy- 
namic self-similarity (DSS). Mathematically, for the scalar decay problem, it amounts to the 
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assumption that 0(r, t) can be reduced to a function of a single variable $(p) through a suitable 
choice of the length and scalar concentration scales i = i{t) and ^ = i9(i): 

e(r,i)=i?2(t)$(^_L_y (2.36) 

The function ^{p) is called the scaling function. There is an equivalent spectral space version 
in which one assumes that 

Eg{k, t) = ^^{t)£{t)F{ke{t)) (2.37) 

for a spectral scaling function F(k). The scale t?(t) can be taken to be the rms scalar intensity 
^(^) '•— [(^^(0)]^ — [0(Oi^)]^ • The length-scale can be taken to be, for example, the scalar 
integral length LQ{t). Within the validity of the hypothesis, all other relevant lengths are either 
0,00, or differ merely by a constant factor. It is natural to take !!.{€) oc Lg{t), because we want 
to consider a limit in which all dissipative length-scales go to zero and we also want to capture 
the energetics of the decay. 

Clearly, the justification of the DSS hypothesis is more difficult than that of PLE. Math- 
ematically, it provides a natural simplifying assumption, but truth and simplicity need not 
coincide. Our main purpose here is to investigate the validity of DSS in the specific context of 
the Kraichnan model. Although the methods we use are very specific to the model, they allow 
us to draw some conclusions that reasonably apply to other problems. The present work should 
also provide a testbed for general frameworks of understanding DSS. For example, it is possible 
that dynamic renormalization group methods can provide a more universal foundation |34|. 



Note: Hereafter energy spectra, integral lengths, etc. will refer to the scalar field only and not 
to the velocity field. Hence we drop the subscript 6 without any possibility of confusion. 
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3 Self-Similar Decay and Its Breakdown 

(3. 1) Derivation & Solution of the Scaling Equation 

In this section we shall completely identify and classify all of the self-similar decay solutions in 



the Kraichnan model. Following the discussion in section 2.2, we look for solutions to (1.9) in 
the form 



As there, we choose 



Q{r,t)=d\t)^[j^^]. (3.1) 



and thus 



r{t) = @{<d,t) (3.2) 



$(0) = 1 (3.3) 



by definition. Subsituting the Ansatz (|3.lD into (1.9) gives 



2^(*) U/„^ , / -^W K^^..^_ 1 9 f^d+c-i9'^ 



Hp) + r. rZu. Pinip) = ^2^in p'^'-'in ■ (3-4) 



yDi^t)LC-2{t) J '■'^' yDiU~^{t)j' dp'-'^' p'^-^dpy dp 

The only way that this can hold with $(/?) a function independent of time t is if 

2i?(t) 



Z)i^(t)LC~2(i) 
and 



-a (3.5) 



^(^) -/? (3.6) 



DiU~\t) 

for some constants a, (3. Indeed, the Wronskian of the two functions ^{p) and p^'{p) in terms 
of the logarithmic variable ^ = Inp is $^(.^)-47 ln<I>(^), so that they are linearly independent 



'd£,^ 



on any interval where $(/>) 7^ and is not a pure power-law. In order that -^{t) < 0, we 
must have a > (which explains our choice of sign in ( |3.5| )). We see that the second equation 
(|3.6|) is identical with ([l.T]) postulated for the scalar integral length, up to the factor of j3. The 



solution is L{t) = [Lq + [3^Di{t — to)] • We may always choose /3 to be unity by a suitable 
choice of lengthscale in the Ansatz. In other words, there are many lengths growing according 
to the equation ( |3.6| ) with some (3 (e.g. the integral length), but we choose the one obeying 
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the equation with (3 = 1. This amounts to a rescaUng /3 —>/?' = 1, L ^ L' = j3^^'^L. This 
transformation takes also a — > a' = a//3. We shall always assume hereafter that this rescaling 
has been done and simply write P' = 1,L' = L,a' = a. Then we arrive at the equation 

The first coefficient equation (p.5|) may then be written using the second (p.6D as 



2^w m 

whose solution is 

T?2(t)=A.L-"(t) (3.9) 

for some constant A with dimensions [0 • L"] . It is clear that any value of the constant A is 
possible, because the equation (|l.9| ) is linear and homogeneous. However, for any solution we 
are always free to chose units of the scalar field (e.g. temperature scale) so that A = \. Then 
( |3.9| ) is identical with ( |1.8D obtained from the PLE hypothesis, if we take a = p. In fact, we can 
verify this directly from equation (p. 7]), if we substitute the asymptotic formula ^{p) ~ Cp~^ 
for p 3> 1. We obtain 

(a-p)p-P~p(d + 7-p)p-(P+^) (3.10) 

for large p. Since 7 > 0, the righthand side is asymptotically negligible compared with the left, 
and we see that a = p. Of course, we do not mean to imply that PLE is necessarily true for 
self-similar solutions, either spatially or spectrally. What is shown above is only that, i/ spatial 
PLE holds, then we have the identification a = p. We may now write out the equation for the 
scaling function in final form, as 

/$"(p) + [(d + C - 1)/"^ + pWip) + a^{p) = 0. (3.11) 

This equation can be solved in terms of confluent hyper geometric functions by means of the 



substitution x = — — . In fact, making this substitution into (3.11) gives: 



dx"^ 



d 
X 

7 



j: $ = 0. 3.12 

ox 7 
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This is a second-order differential equation with coefficients hnear in the variable x. Any 
such equation can be solved in terms of confluent hypergeometric functions. In fact, ( |3.12| ) is 
Kummer's equation with a = a/7, c = d/^ ([^, Ch.6). This equation has two solutions which 
are traditionally denoted $(a, c; x) and ^(a, c; x). The former is an entire function of x defined 

by the power series 

^ / s , ax a(a + 1) x^ , , 

*(a,.;x) = l + -- + -i^- + ... (3.13) 

The second can be defined as 

^(a, c; x) = -Iil^l^$(a, c; x) + ^^^'^K ^-^'^ia - c + 1, 2 - c; x) (3.14) 

l [a — c + 1) l {a) 



when c is not an integer, and otherwise by a limit of this expression (p5[, section 6.5). We see 
that for c > 1 and c non- integer, the second solution has a power-law divergence at x = 0, 
while for integer c the divergence is logarithmic (p5|], section 6.7.1). The Wronskian of the two 



solutions is — w^, (|25], section 6.7) which are thus independent for a,c > 0. These conditions 
are satisfied in our problem since c = - > 1 and a = — > 0. Thus, we may summarize the 



main conclusion of this section as follows: The unique solution of the scaling equation (3.11) 
satisfying the boundary condition <1>(0) = 1 is 

m = ^(-A-^). (S.15) 

V7 7 7 / 
Then, Q{r,t) := d'^{t)^{r /L{t)) with 'd{t),L{t) solutions of ^^ and pT^ j /or /3 = 1 [] is an 
exact solution of \1-^ - Hence, we have identified all possible scaling solutions in the Kraichnan 
model. As an interesting historical note, we observe that von Karman and Howarth already 
reduced their scaling equation, with neglect of triple correlations, to Whittaker's form of the 
confluent hypergeometric equation |j|] . Kraichnan found the energy spectrum of the white-noise 
model to obey Kummer's equation in the Batchelor limit C ^ 2, 7 ^ [16|. Also, in the steady- 



state with random forcing and molecular diffusion, the solution of the 2-point correlation was 
given as a Kummer function by Chertkov et al. in [l^]. 



^If we had used the original length-scale L{t) satisfying (3.6) in defining the scaling solution, rather that L' {t) 

a d. _j3_ 



with j3 set equal to 1, then the result would have been instead $ ' -^ i- —S-rTi 
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(3.2) Asymptotic Behaviors & Permanence of Large Eddies 

The behavior of ^(p) for small and large p can be obtained from the known asymptotics of the 

Kummer function. For small p it follows from ( p. 13 ) that 



$(/,) = l--^p^ + 0(p27). (3.16) 

This is the scaling that we expect for a dissipative solution. Because the Kummer function is 
entire, the only singularity of the scaling function is at p = 0, due to the fractional power 7. 
Hence, we may apply Theorem 4 of Ch.IX.6 of |31] to obtain, for < 7 < 2, that 

$(.) ^ ^ . ^^r f 1±^') r ('^ sin r-1) k^^^^^) (3.17) 



for K —f 00. Hence, from ( 3.16| ) the spectral scaling function goes as 



for K ^ 1. Of course, this is the spectral law one would expect in the inertial-convective range 
of the white-noise model. In fact, going back to dimensionful quantities, we get 

E{k, t) ~ -aC(7, d)t?2(t)L-T(t)K-(^+^) (3.19) 

with C{'y,d) := 2^ \^i/d+'2\ ^i"- (^) ^^^ kL{t) ^ 1. However, referring to the equation ( p. 5] ) 
and using the definition of the scalar dissipation x{t) = ~-dt (i^^(*))' o^^ finds that 

E{k, t) ~ C(7, d)^A;-(i+T) (3.20) 

for kL{t) ^ 1. This is identical to the result that holds in the forced steady state, with the 
same (universal) value of the constant C{'~f,d). In fact, the corresponding spatial result is 

e(r,t)~^2(^)__^(tL_,7 (3.21) 

for r <^ L{t), which coincides exactly with the result in equation (1.19b) of [O. 
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The large p behavior of <&(/)) is obtained from the asymptotics of the Kummer function for 
large negative arguments (|^, section 6.13.1): 

Tic) , 



$(a, c;x) 



r(c-a) 



[1 + 0(1x1-1)] 



(3.22) 



as Rex — > — oo. Thus, we see that 



$(p) ^ ^"/7 






7 



for p » 1. We have set 



v := a — d 



(3.23) 



(3.24) 



so that Q = d+v. We see that ^{p) has the power-law form presumed in the PLE hypothesis, at 
least when v ^ 7^, ^ = 0, 1, 2, .... Otherwise, the coefficient of the above asymptotic expression 
vanishes and it no longer gives the leading behavior at large p. We shall examine that case in 
detail below. However, for all but such exceptional a values the spatial PLE in fact holds. This 



can be seen by returning to the dimensionful variables and using (3^), which gives 

r 



Q{r^t) ~ 7 



a/7_ 



T/ -a 

r 



(3.25) 



for r ^ L{t). This verifies the spatial PLE, for the non-exceptional values a ^ d + ji, i = 
0, 1, 2, ..., since the coefficient of the asymptotic power is explicitly independent of the time. 

Let us now consider the exceptional cases, a = d + ^i, £ = 0, 1, 2, .... It turns out that these 
are given in terms of elementary functions. In fact, for -£ = 0, 1, 2, ... 

1 d' r 



^{c + i,c;-x) 



{c)ix'^ ^ dx^ - 



^c+i-l^-x 



(3.26) 



where (c)^ = c{c+ 1) • • • (c + ^ — 1) and L^ (x) is the generalized Laguerre polynomial of degree 



£. The first line follows from [^], 6.9.2(36) and the Kummer transformation, [^5[, 6.3(7). The 
second line follows either from the Rodriguez formula for the Laguerre polynomial or from p5|. 
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6.4(11) with a = c. Using (p.26|) we see that 



\7 7 7 





(pL^ 


e T . 


\ 7 y 




KfJe 







(3.27) 



ExpHcitly, the first few functions are, for i = 



iy = 0. 



e 1 



for ^ = 1, i^ = 7 



and for ^ = 2, z/ = 27, 



'J>i(p) 



-^ 



_£Z 

e T , 



$2(p) 



1 



2p'> 



+ 



P 



,27 



e T . 



(3.28) 



(3.29) 



(3.30) 



d d{d + j)_ 

We see that all of the scaling functions in the exceptional cases a = d + ji, i = 0,1, 2, ... do 
not behave as power-laws at large p, but, instead, have a stretched-exponential decay. Thus, 
spatial PLE cannot hold. 

The case i = ly = when a = d corresponds to a finite, nonvanishing Corrsin invariant K. 
In fact, the diniensionless version 



K 



^d-i 



00 p7 



can be calculated explicitly by substituting t = p'^ /^ to yield a Gamma integral: 



(3.31) 



K 



^—1 y ( — 

7^ W 

2d-%d/2 r (I) ■ 



(3.32) 



Thus, the spectral scaling function has the expected equipartition power-law at low-wavenumbers: 



F{k) ~ -uJd-iKK' 



d-i _ 7^ W 



.d-l 



(3.33) 



2^-^[r(f)]^ 

Although spatial PLE does not hold, the spectral PLE is valid. Indeed, the dimensionful Corrsin 



integral is 



K{t) = K^\t)L'^{t) = K 



(3.34) 



22 



using ( |3.9D for a = d (with our convention that ^ = 1). The integral is exphcitly independent 
of time, in agreement with the general theorem on invariance. As an aside, we remark that 
the stretched exponential solution was noted by Majda and Kramer [18|. In fact, it was found 



for 7 = 2/3 already by Batchelor in 1952 in a slightly different context [35|. He considered 



L9|) as an equation for the particle pair-separation distribution, in which case it is just the 



scale-dependent diffusion equation proposed by Richardson [36]. However, along with the new 
physical interpretation, there are also different mathematical requirements on the solutions than 
those we impose. To represent a probability distribution function only positive solutions with 
unit integral are admissable. 



Let us investigate the validity of the spectral PLE more generally. It is clear from (3.23) 
that the spectral scaling function F{k) is not C°° at k = 0. In fact, the Fourier transform of 
the radially symmetric function ^{p) is given by 

^('^) = (2^)d/2^(.-2)/2 I P"Mp)Ji,^2)/2{^P) dp (3.35) 



where Jp{z) is the Bessel function |3^. Then, by using the Frobenius series expansion of the 
Bessel function ([^], 7. 2. 1(2)), one obtains the formal Taylor series of <I>(k) in powers of k^: 

oo 

'^{k)^^Bjk^^, (3.36) 

with 

Br-=— -^^7ZP\ p''Hp)p''-'dp. (3.37) 

Because of the power-law decay with exponent —a in ( p. 23 ), one expects that Bj diverges for 



2j ^ i^ 7^ 7^, i = 0,1,2,.... This will be verified below. Hence, derivatives ^^^'{0) of order 
r > iy/2 do not exist. However, it seems reasonable to assume that k = is the only singularity 
of *I*(k), which is C°° elsewhere. In that case. Theorem 4 of Ch.IX.6 of [^l| can be employed 
to infer that the following asymptotic expansion holds for v ^ 2^,7^, j,i = 0,1,2,... and for 
K<C 1 

m 

8(k) ~ Yl Bj^^^ + ^i")'^'' (3-38) 

j=0 
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with 

B{u) = -L- , , ^^' . (3.39) 

^ ' 7r'^/22° r f-i^) r (^) 

In the case u = 2j, z/ ^ ^i, j,i = 0,1,2,..., the second term after the sum is changed to 
B^{u)Ky log K, modified by a logarithm. The summation itself represents the Taylor polynomial 
of degree [^] of the contribution to <I>(k) analytic at k = 0. Its coefficients Bj may be obtained 
from the formula ( |3.37| ) using with j3 = 2j the integral 



p^^i^A--'-] /-Up 



V7 7 7 



r(-)r(^ 



(3.40) 
iboo for P > 1^, 



valid for the nonexceptional values i' ^ ^i, i = 0,1,2,.... The sign it in the second case of 
( 3.40| ) is given by sgnF (— -), the same as the power-law tail in ( |3.23| ). For the exceptional 



values u = 7^, ^ = 0, 1, 2, ... the first line in ( 3.40| ) is valid in both cases f3 < u and j3 > v, with 



the convention that T (-f + 1) /T (-f j = (^-^ + ^ - ij . . . (^-^ + ij (^-^j = (^-^j^. The 
proof of this integral formula proceeds by making the change of variables t = p^ /^ and using 
the Laplace transform of the Kummer function 



00 



-suh-i^i^. „. ^N,, _ r(6)r(c)r(a - 6) _„_. 



e-'H''-'^{a,c;-t)dt = '\ '' ' ,, (1 + sf-^-'F {c- a,l - a;b - a + l;-s) 
r(a)r(c-6) 

+ ^^^l^^^~°^ g°-'F(a,l + a-c;a-6+l;-g) (3.41) 
1 (c — aj 

which is given in terms of the hypergeometric function F and valid for Res > 0. This formula 
follows from |25],6.10(5) and 2.10(4). Notice that for the exceptional cases a = c + i, £ = 



0,1,2,... the second term in ( |3.41| ) vanishes. The limit s J, may be obtained recalling that 
F{0) = 1, with the result given in ( |3.40| ). We see that in general for v > 0, there are nonzero Bj 
terms from the analytic contribution which dominate the Fourier transform at small k. However, 
it is interesting to note that the lowest coefficient Bq, which coincides with the dimensionless 
Corrsin integral K, always vanishes for z^ > 0. Thus, the Fourier transform $(k) always vanishes 
at K = for z^ > and the second singular term proportional to k'^ is formally the leading-order 
contribution for all i^ < 2. 
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We may now return to the question of spectral PLE. The asymptotic low-wavenumber 



expansion for the spectral scaling function, corresponding to (3.38), is 



F{k) ~ Yl AjK'^+^^-^ + ^(i/)k"-i (3.42) 

j=0 



valid for k <C 1, with Aj := -^uJd-iBj and 



A{u) := -^ y^ • ,,y^, , . (3.43) 

2- ^f_.^ r f r f 



If we return to dimensionful variables and use again (|3.9| ), we obtain 

m 

i=o 
for kL{t) <C 1. Spectral PLE appears to hold, in the sense that the coefficient A{v) of the 
singular term is explicitly time-independent, for all but possibly the exceptional values v = 
2j,^i, j,i = 0,1,2,.... Of course, the singular term is the leading one for ly < 2 and then 
spectral PLE appears to hold in the standard sense. 

However, this is wrong for an important reason. While A^v) > for — d < i^ < 7, it becomes 
negative at rather higher values of v. A{v, d) < for 7 < z^ < 2 min{l, 7}! Thus, realizability of 
the scaling solution is violated in this case. This is a crucial issue which we have neglected up 
until now. In fact, only solutions with positive spectra over the whole wavenumber range are 
physically admissable. Thus, we have reached one of the important conclusions of this work: 
No self-similar decay is possible in the Kraichnan model with 7 < i^ < 2min{l,7}, since in 
that range the scaling solution has a negative spectrum at low wavenumbers. We have proved 
this subject to a single assumption, that the origin is the only singular point for the Fourier 
transform of the scaling function. The result is verified by explicit computations for the special 
case 7 = 1 in an Appendix. In fact, we shall prove below that self-similar decay occurs for no 
exponent i^ > 7. The physics of this phenomenon will be discussed in section 3.4 below. 
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(3. 3) Realizability of the Scaling Solutions 

We must now examine more closely the issue of realizability. As already shown by Kraichnan Q , 

the 2-point correlation of the statistical problem (1.3) must satisfy the closed partial differential 



equation ( |1.9| ). However, there can be solutions of the PDE which do not correspond to any 
solution of the statistical problem. The necessary and sufficient condition for a solution of 
( |1.9| ) to be realized as a solution of the statistical problem ( |1.3| ) is that it be positive-definite. 
Necessity is obvious. To see sufficiency, take for any positive-definite initial data 0(r, 0) the 
Gaussian measure ^o over scalar fields ^o(r) with zero mean and with the given positive-definite 
function as its 2-point covariance. Then, the solution 0(r, t) of the PDE (|1.9| ) with the specified 
initial datum 0(r, 0) will be the same, by Kraichnan's result, as the 2-point correlation of the 



statistical problem posed by (1.3) with random initial data ^o(r) distributed according to fiQ. 
It is a corollary of this remark that the PDE (|1.9| ) is positive-definiteness preserving, that is, 
positive-definite solutions result from positive-definite initial data. The upshot is that not every 
scaling solution of ( |1.9D that we found in Section 3.1 can necessarily be realized as a solution 
of the Kraichnan model ( |1.3D . This will only be true if it is positive-definite, and this question 
must now be addressed. 

First, we give a general proof of realizability when 0<a<(ior— (i<z^<0. We must show 
that ^(p) is positive-definite as a function on d-dimensional Euclidean space. This is proved in 
two steps. We observe first that the v = function <&(-,-;— ^j =e ~ is positive-definite 
for < 7 < 2. In fact, for < 7 < 2 the functions e t are the characteristic functions (Fourier 
transforms) of positive probability densities, the spherically symmetric stable distributions of 
parameter 7. General multivariate stable laws were first investigated by Levy |3g] and Feldheim 



|39|; for an introduction to their basic theory, see, for example, the monograph of Zolotarev 
1 40|, Section 1.6. The positive-definiteness of the characteristic functions e ^ for the spherically 
symmetric stable distributions can be obtained by an easy modification of the proof of Bochner 
for the 1-dimensional case pl|. However, with this result, we may then use a standard integral 
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representation for the Kuninier functions 

Ha,c;x) = —^ ^ [\-^u--\l-ur—'dn (3.45) 

T{a)T[c - a) J^ 

valid for Rec > Rea > 0. See [25|, 6.5(1). From this we see that 



$(—,—; = — / \ / \- / e ^ u^ (1 — u) T an, (3.46) 



7 7 7/ ri^r-^Wo 



when —d < u <^. Hence, <!>(—,-;——) is a convex combination of positive-definite functions, 
in fact, averaged over a Beta distribution. Thus, this scaling function is always positive-definite 
for -d < z^ < 0. 

We now show that the scaling functions are positive-definite — and thus realizability holds — 
for the larger range —d < i^ < 7 or < a < (i-l-7. As we have seen above, the low-wavenumber 
spectrum remains positive in the range —d < z^ < 7. Furthermore, we prove now that, if the 
scaling function for z^ = 7 is positive definite, then so are all the functions for —d < z^ < 7. 
This follows from the identity R 

^^' + ^' "' "^^ = r(c + ^')ra-^) r ^^' ^ ^' "' ^'^^"'^'"'(^ " ^)~''^^ (^-^^^ 

for 1 > Re^ > —Rec. This is most easily proved by expanding both sides in a power-series in x 
and comparing the coefficients. One may also give a proof based upon [p3], 6.4(12) and 6.5(1). 
Thus, 

^{d + V d p^\ ^VV') f^^fd + J d p^'\ i±JL-X, ,^^ , . . ^ 

$ ,-;-—]= —^ \—^ ^ $ -,-;-u—]u'y {I - u) '' du, (3.48) 



7 7 7/ r(^^±^)r(3^^)^o V7 7 7 

when —d < 1/ < 7. Hence, if the scaling function for z^ = 7, <^ ( — -, -;——), is positive- 
definite, then, as weighted integrals of it with respect to a Beta distribution, so are the scaling 
functions for —d<v<^. 



^Although it is not hard to prove, we did not find this formula in standard treatises on Kummer functions. 
We would be grateful for any reference. 
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It therefore becomes important to answer whether $i(/j) = ^ i- + 1, -; —^ I is positive- 
definite or not. We shall show that this question is related to properties of the equipartition 

d d. p''' 

- - J . ill XCH^U, lU lO IIWL IK 

d 



solution $o(/o) = ^ ( -, ^; — T7 ) • In fact, it is not hard to show by direct calculation that 



d-<^i{p) 



d + p- 
dp 



Mp)- (3.49) 



This is actually a particular example of a standard relation between Kummer functions, |25|, 
6.4(11), transformed by the substitution x = —p'^/j. What is useful here is that the righthand 
side of ( p. 491 ) above is also just minus the lefthand side of ( |3.7| ), in the equipartition case a = d. 
However, the expression in ( |3.7D arose by differentiating the scaling Ansatz with respect to time 
and dividing by DiL~'^{t). Thus, we infer that 

where &o{r,t) = ^^{t)^o{r/L{t)). This relation may be Fourier-transformed easily, giving 



If we employ the relation QQ{k,t) = ^Q{kL{t)), which follows using ( p^ ) for a = d, and ( p^ 
for the time-derivative of L{t), then (|3.51) yields the final result 



$i(k) = -^^^{^)- (3.52) 

d dn 

From this we can see that realizability holds for the threshold case ^ = l,i/ = 7 when the 
Fourier transform ^q{k) is a monotone nonincreasing function of spectral radial coordinate n. 
To complete the proof, we must verify this property. In fact, the monotone nonincreasing 
of the density with respect to the radial cooordinate is equivalent — for spherically symmetric 
functions — to the property of unimodality of a multivariate density, as it has been defined by 
Olshen and Savage ^^. For a general introduction to the subject of unimodality and to its 
proof for stable probability laws in particular, see pO|], Section 2.7, and p^]. The subject has 
a rather colorful history, involving a series of published false proofs and claims by eminent 
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mathematicians (including Kolmogorov), which is summarized in those works. The first proof 
of unimodahty of the symmetric, one-dimensional stable distributions was given by Wintner 
in 1936 |44]. It was then widely conjectured that all one-dimensional stable distributions are 



unimodal, but it took over forty years until a correct proof was found in 1978 by Yamazato 
|45|. The proof of unimodality of spherically symmetric stable distribution functions in multi- 



dimensions was given about the same time, by S. J. Wolfe [^]. This is exactly the property we 
need to guarantee realizability for the threshold case u = j, and thence, by equation ( |3.48| ), for 
all z^ in the range —d < v < j. 

On the other hand, it is reasonable to conjecture that no scaling solutions with u > j are 
realizable. We have already established that realizability fails when 7 < z^ < 2min{l,7}, by 
showing that the low-wavenumber spectrum becomes negative. The same argument does not 
work for all 1/ > 7, because (i) the coefficient of the singular contribution to the spectrum in 
(|3.42| ) oscillates in sign as u is increased and (ii) the singular term is not the leading-order 
term at low-wavenumbers for zv > 2. However, we shall now prove that the scaling solutions for 
1/ > 7 are indeed non-realizable. The proof is based upon the fact that, for z^ > 7, two integrals 
vanish, namely: 

f ^p)(fp= f pf^{p)d'^p = 0. (3.53) 

This is a direct consequence of the general formula ( |3.40| ) for the case /3 = 7. We shall now 
show (following a suggestion of D. Thomson [47, ^]) that these two conditions are equivalent 
to 

^-(d+7)$(^)^rf^ = 0, (3.54) 



/ 



when realizability (or nonnegativity of ^t pointwise) is assumed. To prove ( 3.54| ) we make use 



of the following generating functional for the moment-integrals in (|3.53D 



1 / -t^. 



G{t) ■■= T^TTR / e-'-^p)d^p. (3.55) 

(2vr)'* J^d 



Then it is easy to see that (3.53) is equivalent to G'(O) = 6(0) = or, as well, limt^o ~i = 0- 
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Next we make use of Parseval's theorem to rewrite ( |3.55| ) as 

where $o(^) is the Fourier transform of the scahng function for i^ = 0, i.e. the density of the 
multidimensional Levy stable distribution with parameter 7. This step is justified because, 
clearly, <I>o ^ L'^i and because, for general i^ > 0, boundedness and the large-p decay in ( 3.23 ) 
imply that as well <I> G L^. Next we observe from ( 3.17] ) that $o('^) ~ c • k"'^'^^'^' as k ^ 00, for 
some positive constant c > 0. Thus, 

lim ,J,, $0 (4r) = c • '^'^'^^^^ ■ (3.57) 

We see finally, by Fatou's lemma, that 

= liminf^ = liminf / ,_\, $0 f-^") $(k) d'^/^ 
> I liminf ,,^\, ^o(^)^{K)d'^K 
= c/ K-('='+^)$(K)d^K>0, (3.58) 



where the last inequality holds by the assumed nonnegativity of <1>(k). Thus, the identity (3.54) 
is established. However, it follows immediately then that <5(k) = for a.e. k, which is a clear 
contradiction. Thus, the assumption that <1*(k) > for all k cannot be correct. 

We may summarize the conclusions of this section as follows: The scaling functions ^{p) = 
<!>(—,-;——) are positive-definite, radially- symmetric functions on M for all a in the interval 
(0, d + 7] and for no a in the range {d + 7, cxd). In particular, the scaling solutions 0(r, t) = 
•d"^ {t)^{r / L{t)) for a < d + 7 are all realized as solutions of the statistical problem posed by 
the equation ( [j.^ with the random initial data fiQ. It is interesting to note that, for cases 
z^ > 7 when the low-wavenumber spectrum is positive, there must be at least one negative 
spectral interval away from the origin. This is explicitly verified in the Appendix for 7 = 1, in 
which case the spectral scaling functions F{k) can be calculated in a closed form (related to the 
multidimensional Cauchy distribution) for all of the exceptional values 1^ = 1, i = 0,1,2, .... 
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(3.4) A. Physical Explanation of the Results 

We have now shown there are reahzable self-similar solutions for —d < v <"f but not for z^ > 7. 
Although v = J corresponds to a realizable DSS solution, something a bit strange must occur 
in that case. For example, PLE cannot hold either spatially or spectrally. If there is a power- 
law in Ei{k,t) at low-wavenumber, it must be distinct from the power k'^^"'^'^ which is naively 
expected. Indeed, if that naive power-law occurred, then it would imply a corresponding spatial 
decay r~^'^~^'^\ since 7 lies in the range < 7 < 2 where this deduction is correct. However, 
we know that 0i(r, t) decays faster than any inverse power of r. Hence, if any power-law at all 
occurs in the spectrum, it must be different from the naive one and, in fact, it must correspond 
to V an even, positive integer. We show now that there is a spectral power-law k*^^^ (naively 
corresponding to 1/ = 2) and representing scalar backtransfer. 
More precisely, we show that 



$i(At) 



d+i -p ( d+2 
1 '' V ^ / ^2 



(3.59) 



^(4^)d/2 r (d±2^ 

for K <^1. To prove this, we use ( |3.51| ) again, but calculate now the time-derivative using the 
transfer equation discussed in Section 2.1: 



dtQoiht) 



rCiKj 



AAj(q) Go(A;,t)-eo(|k-q|,t) 



(3.60) 



where -Dij(q) = K{q)Pij{c\). We can take the limit as /c ^ in this expression, with the result 
that 



dtQo{k,t) 



d-1 



d 



k' / (i^qK(g) eo(0,t)-Go(g,t) 



(3.61) 



where the factor {d — l)/d comes from the angular average of Pjj(q). Fourier inversion gives 

J A K{q) [Go(0, t) - Qo{q, t)] = -^ I d'^v [K(0) - K{r)]Qo{r, t). (3.62) 

In this form, the limit as /cq ^ is easy to evaluate, with the result 

lim f d'^Y [K{Q) - K(r)]eo(r, t) = Ki [ d^r r^eo(r, t), (3.63) 

fco^O J J 
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for Ki = ( ^^ j Di. By substituting the scaling Ansatz and using again (|3.9|) , one finds 

//•CXD 
d-^r r^eo(r, t) = L«(t) • 00^-1 j p'^+^-^^^ip) dp. (3.64) 

The right side of ( p.64 ) can be reduced to a Gamma integral by the subsitution t = p'^ /^, giving 

/+^-i<5o(p) dp = t'^-^F ('^-±^\ . (3.65) 



Putting together ( 3.61[ )-( pl65 ) gives finally 

/in 11^ A ^1 fd-l\ (d + Q 

dt@o{k,t) ~ 



(27r)'>' V d J \d-l 



t^d-i7 ^ 



^-^rf^^)LC(t)(./L(t))^ 



^-^^ + 2^ W(4.)'^/2r^^ + ' 



(3.66) 



7 7/ V 2 

Substituting this into ( 3.5l|) yields ( 3.59 ), as claimed. 

The result may be easily understood in terms of our general result ( [j.38 ) for the low- 
wavenumber asymptotics. The coefficient B{u) of the singular contribution vanishes at 1/ = 7, 
because of the Gamma function Fl— M in the denominator of ( 3.39|) . Thus, the leading 
term ought to be Bik?. If we substitute /3 = 2,z/ = 7 in ( p. 401 ), the integral is found to be 



Ir ( d+2\ 



d+2 



7 



. This yields Bi 



d+( p/ d+2^ 
7 ^ I 7 y 

d+2' 



d(4Tr)d/2 p^ 



in ( p. 371 ), in exact agreement with ( 3.59| ) 



In terms of the spectral scaling function, the result for our DSS solution is 



F,{k) 



d+l-i -p ( d±2 
77 V 7 



^d+1 



(3.67) 



2^+1 [r(^)]^ 

for K <C 1. Although self-similar, spectral PLE must be violated in this solution. In fact, 
when the naive power is replaced by any other, the DSS solution will automatically develop a 



time-dependent low-wavenumber coefficient, in contrast to (3.44). To see this, we may return 
to dimensionful variables in ( p. 671) , using now ( |3.9|) with a = d + j. The result is 



Ei{k,t) ^A{j,d)L^{t)k'^+^ 



(3.68) 



-1 r 



d+2 ' 



for kL{t) <C 1, with A{'y,d) := "^^^^^ }d+2\]^ - Because of the "leftover" factor of L'^{t), 
the coefficient of the asymptotic k'^^^ power-law is now explicitly time-dependent. In fact, it 
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exhibits power-law growth L''{t) ~ [^Di{t — to)] for long times t — to ^ L'^{tQ)/Di. This 
increase is consistent with the interpretation of the low-wavenuniber spectrum as arising from 
scalar backtransfer. 

Since there are no realizable self-similar solutions for ly > j, one is led to the following 
question: how shall an initial scalar spectrum Eik^to) ~ Ak'^°~^ at small k with oq = d + uq 
decay asymptotically at long times? According to the traditional view, when vq > 2, then 
the decay shall be asymptotically self-similar at long times described by the spectrum Ei{k,t) 
for the threshold case v = j above. We see no reason to doubt the validity of this view. 
On the other hand, when —d < uq < 2, the traditional view states that the decay will be 
described at long times by the self-similar solution with v = uq. This is perfectly consistent 
when —d < z^q < 7 and, again, we see no reason to doubt the traditional picture of the decay. 
However, when 7 < i/q < 2, then there is no realizable self-similar solution with v = uqI Hence, 
the traditional view must be wrong for 7 < i/q < 2. 

We cannot from our present analysis say what happens in the case 7 < 1^0 < 2, because it 
lies outside consideration of the self-similar solutions themselves. We shall treat this dynamical 
problem in Section 4. However, we shall now show that the scenario proposed in Q for the 
Burgers decay gives a completely consistent account of the known facts also in the Kraichnan 
model. The picture that is proposed is of a two-scale decay. In addition to the integral length 
L(t), there is another length-scale L^[t) ^ L[t), which separates an inner solution Ein{k,t) 
for kL^{t) ^ 1 and an outer solution Eout{k,t) for kL^:[t) <C 1. The inner solution is just 
the self-similar decay solution Ei{k,t) for v = j. The outer solution is the same as the initial 
spectrum Eout{k,t) ~ ylA;""~^ for kL^{t) <^ 1. Hence, spectral PLE holds, but only in the outer 
region. Now by matching the inner and outer solutions, one can find the crossover length-scale 
L*(t), or, equivalently, its associated wavenumber A;^=(t), as 

A[K{t)Y+'''''^ ~ D{t - tQf'^%{t)Y+\ (3.69) 
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or 



L^{t)r^{t-to)(^-'^-o)-t. (3.70) 

Since L{t) ~ (t — to)^ ) the inequality L^,{t) ^ L{t) necessary for validity of this picture only 
holds if C/(2 — t'o) > 1 or if 2 > i/Q > 7- It is a highly nontrivial test of consistency that the 
critical value uq = 7, above which the ratio R{t) := L^{t)/L{t) grows, coincides exactly with the 
value z/ = 7, above which there is no realizable self-similar decay solution. This gives us some 
confidence in the correctness of the picture proposed. If this picture is correct, then at very 
long times DSS is restored and the decay is described by the inner solution spectrum Ei{k,t), 
since -L*(t) — s- +00 in units of L(t). As fo ~^ 2—, the growth rate of L^(t) becomes infinitely 
fast and the outer solution region disappears, in agreement with the traditional view. Thus, 
PLE and DSS both hold for 7 < i/q < 2, but the decay is not what one would naively expect 
for DSS+ PLE, because the outer range where PLE holds is not part of the inner-range DSS 
solution. 

A similar picture may be developed in physical space, but, as noted in O, the separation 
of the length-scales is not as sharp. Now one would expect 0m(r, t) = @i{r,t) for r <C L^(t) 
and @outir,t) ~ AV-('^+''") for r » L^t). Note that for r ^ L{t), 

eUr,t) ~ [L(t)]-('^+^) » Qout{r,t) ~ [L(t)]-('^+'^o). (3.71) 

The crossover occurs at a larger length-scale L^{t), which is found by the matching condition 

[L(t)]-('^+^) exp[-(l/7)(L,(t)/L(t))^] ~ [L,(t)]-('^+^o). (3.72) 

It is easy to see that this implies a solution for the ratio i?*(t) = L^{t)/L{t) of the form 

R{t)^[log{t-to)]'/\ (3.73) 

Hence, L^{t) is only larger than L{t) by a logarithmic term in physical space. 
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4 Convergence to Self-similar Solutions 

(4-1) The Long-time Scaling Limit 

In this section, we identify the time-dependent solutions of the Q equation ( |1.9| ) which show 
eventually a self-similar form of decay. That is, we find domains of attraction of the self-similar 
solutions constructed in Section 3. We are interested in observing the solutions on a range 
of length-scales comparable to L{t) and at a level of scalar amplitude comparable to the rms 
fluctuation i?(t). Hence, we consider the rescaled solutions 

e(r,t)=^2(t)r(r/L(t),r(t)), (4.1) 

where | = -f ^, ^ = DiL'^-^ are as in (|3j),(U) and 

r(t):=logL^~log(t-to). (4.2) 



The function F solves: 



dr^ ^ Id 



.-<-|(.,.) 



dT 

+ p—{p,T) + ar{p,T), (4.3) 



5r ' p'^ ^ dp 

where p = r/L{t). This change of variables may be made for any a > 0, and, if we wish to 
make this explicit we shall refer to the above function as '^(a) (P) ''")• Of course, it is only possible 
that a nontrivial scaling limit is obtained as r — > cxd for one exponent a. If a nontrivial limit is 
obtained for a, then, for any other exponent, say, a', it follows that 

oo if a' > a 



lim r(„,)(/9,r) = < 






(4.4) 
if a' < a 



Of course, this means that if the scaling limit Ti^t\{p,T) -^ (resp. oo) for q', then an exponent 
a larger (resp. smaller) than a' is required for a nontrivial limit (if this is possible at all). 

To analyze the equation ( [4. 3D , we make the change of variables x = ^~^p'^ ^ which is the 
same as in Section 3 up to a sign. We obtain 

dV d'^r dV 

— (x, r) = x-^{x, t) + {c + x) — {x, t) + aT{x, r), (4.5) 
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with a = a/7, c = fi/7. Let $(q)(2;) := ^{a,c;—x) be the self-similar solution discussed in 

Section 3, for the same choice of a as employed in Tj-^) . Then $ is a steady state solution of ( |4.5[ ). 

We define A{x,t) via T{x,t) = ^{x)A{x,t). Clearly, to establish convergence T{x,t) -^ <J>(x), 

it is enough to show that A{x,t) — > 1. We find for A that: 

dA^ , d'^A^ , / <^'{x)\ dA^ ^ 

57(^>^) = ^5^(^>-)+^c + ^ + 2x^^(x,r) (4.6) 

The transformation to A(x, t) has removed the explicit a dependence, which is now represented 
only through ^[a)- 

We shall now obtain characterizations of the domains of convergence of the self-similar 
solutions for each a. In this analysis, it is important to distinguish two general classes of initial 
data: those of "rapid decay" for which 

|e(r,0)pe'^(ir < 00 (4.7) 

and those of "slow decay" for which the above integral is infinite. Thus, the initial data showing 
rapid decay belong to L^ with a stretched-exponential weight and must decay at least as fast 
as the weight. This turns out to be a useful formal criterion for "rapid-decay". Initial data 
with power-law decay at large r — which are of particular interest in view of the question of 
validity of spatial PLE — are classified as "slow decay" functions. Convergence results will be 
established for both classes of initial data below. 
(4-2) Initial Data with Rapid Decay 

We analyze first the case of "rapid decay". If we make the scaling with a = d and define 
correspondingly T{x,t) = A{x,t)(^q{x), then ( ^ ) becomes 

— (x, r) = x—-^{x, t) + {c- x) — {x, t) := CoA{x, t) (4.8) 

We used ^o{x) = e~^. The operator Co has as its eigenfunctions the generalized Laguerre 
polynomials L^~ (x) with eigenvalues —£. Indeed, the Laguerre polynomial is characterized as 
the unique solution y = L^~ (x) of the second-order equation 

xy" + (c - x)y' + ey = (4.9) 
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which is regular at the origin. See ||25|, 6.9.2(36) or |^, 5.1, 5.3. In terms of A{x,0), the rapid 
decay criterion ( |4.7| ) becomes 



r-oo 

|2„-a- c-1. 



/■c 

/ |yl(x,0)|^e-^x'=-^dx<oo (4.10) 

Jo 



It is known that the generalized Laguerre polynomials {L'j,~ (x) : £ = 0, 1, 2, ...} are a complete, 
orthogonal set on the interval (0, oo) with the weight Wc{x) = x'^^^e~^. See ^^, Theorem 5.7.1. 
Thus, it follows from the above remarks that the solution A{x,t) of (4.8) with initial datum 



A{x, 0) also satisfies the condition ( [4.1[1| ) and, furthermore, has the expansion 

oo 

A{x, t) = Y, aie-^^ L^f^x) (4.11) 

e=o 

which converges in the L^-sense with weight Wc{x). The expansion coefficients are given in 
terms of the initial datum by 

ai = TTT—^ / Aix,0)L''f\x) e~^x^-^dx. (4.12) 

r(c + ^)yo 



We assume the standard normalization of the Laguerre polynomials, |4£], 5.1.1. 

An important conclusion follows immediately from the fact that the expansion coefficients 
may also be obtained from the solution A{x,t): 

«^e"'" = FT^^ / A{x,t)^-\x) e-^x'-'dx. (4.13) 

i~[C + I) Jo 

If we recall that T{x,t) = A{x,T)e~'' and that e^ = L^(t), then (|4.12D and (|4.13| ) together 
imply that 

/•oo /•CX) 

[L{t)Y^ T{x,T)Ll-\x)x^-^dx= V{x,Q)Ll-\x) x'^-^dx. (4.14) 

JO ^0 

In other words, the lefthand side of this equation is for each £ = 0, 1,2,... an invariant of motion 
of the equation ( [4.5| ) (with a = c). If we return to the unsealed solution, these invariants, with 
an appropriate choice of normalization, take the form 
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for ^ = 0, 1, 2, .... Thus, there is an infinite sequence of integral invariants of the equation ( \1.S( ) 
for Q{r, t). The first such invariant for £ = is nothing but the Corrsin invariant Jq = K. The 
next two are, for i = 1, 



Mt) ■= 



^""W - 1 



e(r,t) d'^r 



and for i = 2, 



Mt) ■= 



i2.(,)_?B^ + 



.27 



Q{r,t) (Ty. 



(4.16) 



(4.17) 



d ' d{d + 'j)_ 

These are really "generalized invariants", because they depend not only upon the solution 

Q{r,t) but also explicitly upon the time t. However, we see that, in the subspace defined 

by the vanishing of the first p invariants, Jq = Ji = • • • = Jp-i = 0, the pth integral Jp oc 

j r'P^Q{r,t) d'^r and is an ordinary invariant. This may also be inferred directly from equation 

(|1.9|) , by using the fact that its righthand side defines an operator homogeneous of degree —7 

and using integration by parts. 

These new invariants play a key role in the problem of the convergence in the scaling limit. 

We see using the relation T{x,t) = A{x,T)e~^, the definition of ^i{x), and the expansion 

formulae (|1|),(|1|) that 

00 

T{x,T) = J2ci[L{t)]-^'^e{x) (4.18) 



e=o 



with convergence of the summation in the L'^-sense on (0, cxo) for the weight Wc{x) = e^x' 



Xr^C—l 



The expansion coefficient 



1 



C£ 



T{x,0)L''f\x)x^-'dx 



(4.19) 



r(c) Jo 

is proportional to the invariant Je, i = 0,1, 2, .... Now suppose that the first p invariants vanish: 
Jo = Ji = ■ ■ ■ = Jp-i = 0. Since the solution $p(x) corresponds to a = d + P7, we see that if 
we scale according to that a, we obtain 



Tix,T) = Y,cp+dmr^''^p+ei^ 



(4.20) 



e=o 



Hence, it follows that lim.r^oo^{x,T) = Cp^p{x). Gathering together the above results, we may 
state the following proposition: Suppose that the initial datum 0(r, 0) lies in the Lp' space with 
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stretched exponential weight e t and that the first p invariants vanish Jq = Ji = ■ ■ ■ = Jp-i = 
but Jp ^ 0. If one scales the solution Q{r,t) with a = d + pj, then 

lim r(„)(x,r) = Cp<l>p{x) (4.21) 

with convergence in the L^ sense on (0, oo) with weight Wc{x) = e^x'^~^. The constant Cp is 
given by the ratio of the pth invariants Jp of the initial datum T{x, 0) and Jp of the equilibrium 
solution ^p{x). It only remains to justify the last claim. In fact, with the normalization of the 
invariants adopted here 

fOO 



1 f°° 
Jp = :=--- / ^p{x)Ll-^{x)x''-'^dx 



= 1. (4.22) 

This completes the proof. 

It is interesting to note that Jq = Ji = but < Jp < oo for some p > 2 is not consistent 
with realizable initial data, due to the nonrealizability of scaling solutions ^p{x) for p > 1. 
Indeed, if p is the least integer p for which Jp ^ 0, then our preceding result implies a fortiori 
that r(rc,r) -^ Cp^p{x) in L^ with respect to the finite measure e~^x'^~^dx and hence, along 
a subsequence of times r^, k = 1,2,..., convergence for a.e. x. If the initial data were realiz- 
able (positive-definite), then, since positive-definiteness is preserved by the dynamics and by 
pointwise limits, the limit Cp^p{x) would be positive-definite as well. However, this contradicts 
our earlier result. Thus, no positive-definite 0(r, 0) can have Jq = Ji = but < Jp < oo for 
some p > 2. Our argument here is rather indirect, using the equation (|1.9D , but the conclusion 
involves no dynamics. In fact, it follows directly by the same argument used in section 3.3 to 
prove nonrealizability of the scaling solutions for i/ > 7 that any initial datum 0(0) G L^ (un- 
weighted) with Jq = Ji = cannot be realizable (positive-definite) . The condition Jq = Ji = 
is precisely equivalent to the condition of vanishing moments, ( p. 531 ), employed there. 
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There is a direct connection (pointed out to us by K. Gaw§dzki) of our generalized invariants 
with the "slow modes" in 2-particle Lagrangian statistics that were discovered in ^^. It is 
a consequence of that work that there are homogeneous moment functions (j)o^p{r) of degree 
f^'o.p) P = 0,1,2,..., whose integrals over any initial 2-point function evolve in time as pure 
degree p polynomials: 

f <Ao,p(r) e(r, t) d^v = Y, Cp,q tP-' I hA^) e(r, 0) d^r, (4.23) 



for some computable constants Cp^q = /'0o,(j(r, 1) (j)o,p{r)d'^r, in the notations of [^]. As these 
constants are manifestly independent of initial data, it is not hard to infer from (4.23) the 



existence of an associated sequence of "generalized invariants". These "slow modes" were 
constructed in |5C] for every angular momentum sector i = 0, 1, 2, ..., but for the rotationally- 
invariant sector 1 = they are particularly simple, given just by the powers (po,p{f) = f^- See 



Appendix A of |5C]. It can then be easily shown that the associated sequence of generalized 
invariants in the i = sector coincides with the sequence Jp, p = 0,1,2, ... we found above. 
(4-3) Initial Data with Slow Decay & A Finite Invariant 

The previous results do not allow us to address the question whether permanence of large eddies 
(PLE) holds in the space domain. For this, we must consider initial data with only power-law 
decay at large distances. Such data with slow decay fall themselves into two broad classes: 
those with one of the invariants Jg, i = 0,1,2, ... finite and those with no finite invariants. Here 
by "finite" we mean both non-zero and non-infinite. It will be shown below that the class of 
initial data with a finite invariant behaves very similarly — as far as the leading-order behavior 
is concerned — to the initial data with rapid decay. 

To study such initial data, a new technique is required. An important observation is that 
(|4.6| ) is the backward equation corresponding to the Fokker-Planck equation on the half-line 
x>0 

^{x,t) = -^ {a{x)P{x,r)) + ^^ (6(x)P(x,r)) (4.24) 
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with drift 



and diffusion 



a(x) = c + X + 2x 



$'(x) 
$(x) 



b{x) = o"^(x) = X. 
Thus, it follows immediately that 

A{x, t) = E[AiX,^r,0)] = / dy A(y, 0)P(y, t\x, 0) 

Jo 

where X^^r is the diffusion starting at x at r = and obeying the Ito equation 

dXr = a{Xr)dT + a{Xr)dWr, 



(4.25) 



(4.26) 



(4.27) 



(4.28) 



and -P(y, t|x, 0) is the transition probability kernel for the process. The formal invariant density 
P{x) of the Fokker-Planck equation is 



P(x) ocx^-^e'<^'^(x). 



(4.29) 



However, this is only a normalizable probability density for the exceptional cases when v = 

7^, I = 0,1,2,..., in which case we refer to the density as Pi{x), i = 0,1,2,.... This can be 

easily understood from the character of the drift term. When <l>(x) oc x'"' at large x, then 

a(x) ~ c — 2a + X for x ^ 1, which is unstable. However, when <l>(x) = p(x)e~^ for some 

polynomial p(x) — as for the exceptional series — then a(x) ~ a — x for x ^> 1, which is stable. |^ 
^These observations give another method to infer the existence of the sequence of generahzed invariants 



Jp, p = 0, 1, 2, ... In fact, using the invariant density in (4.29), equation (4.6) may be rewritten as 

d 






x-e""^ (x)—-(x,t) 
ox 



It then follows by a formal integration by parts argument that 

/"OO 

J{r) ■- / dx a:""'e"yl(a:,T)$2(x) 
Jo 

is an invariant for all a. However, only for the exceptional values a = d + 7p, p = 0, 1, 2, ... do the integrands 

decay rapidly enough to justify this argument. These give the familiar sequence of invariants Jp, p = 0, 1, 2, ... 
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In particular, it is true for the equipartition solution ^o{x) = e ^ that 

a{x) =c-x (4.30) 

and the invariant distribution is given by 

Poix) = ^^x'-'e-^. (4.31) 

For this i = process we can also calculate the transition probability kernel by means of 
expansion in the eigenfunctions of the backward operator Cq in section 4.1: 



Pfe..|..„) = i_^j:iLWi^,_. „,3,) 



r c ^^ i+c-i 



This sum can be evaluated in closed form, using |g5|], 10.12(20), with a = c— l,z = e '^: 



in terms of the modified Bessel function /c_i(z). Cf. |13|, equation (2.15b). 

Using these results, we can establish our first main convergence result of this section: // 
0(0) € L°°P|L^ and the Corrsin invariant of the initial datum is non-zero, then the limit 
exists, 

lim '^(d){x,T) = co^oix) 

T— »00 ^ ' 

uniformly on compacts in x. The constant cq is the ratio of the Corrsin invariants K of the 
initial data r(x, 0) and K of the equipartition scaling solution ^q{x). The condition B(0) E L°° 
is natural, since |G(r, 0)| < 0(0,0) < cxo for any positive-definite initial data with finite energy. 
This theorem includes the case of power-law decay 0(r, 0) ~ Ar~" for r ^ Lq with a > d, 
which guarantees integrability. Because 0(0) G L^, the Corrsin invariant must be finite. The 
theorem remains true even if the Corrsin invariant is zero initially, but in that case it yields a 
trivial (null) scaling limit. 

We prove the result first for bounded ^4(0). This will follow from standard convergence 
results for time-dependent distributions of one-dimensional diffusion processes. For example, 
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T[, Section 23, Theorem 3 implies that for any A{0) € L°° and uniformly on compacts in x, 

hm E [A {X,,r,0)] = [ A{y, 0) P{y) dy (4.34) 

where 

P{x) ex -^ exp I /' ^dtj . (4.35) 

is the invariant measure of the process, if the diffusion is defined with instantaneous reflection 
upon hitting the boundaries. Such boundary conditions must therefore be checked to be sat- 
isfied. We show, in fact, that the process X^^r has zero probability of reaching either of the 
boundaries of the semi-infinite interval, or oo, at any finite time. Thus, it is a special case of 
reflected b.c, with no reflection ever required, and the Theorem 3 of pll] applies. 

We will treat the boundary conditions in a generality that will permit to discuss later cases 
as well. According to [51|, Section 21, Theorem 1, for any b G (0,oo), if 



L- = / exp i - / -^rAdt \ dx = -Foo, (4.36) 

Jo [ Jb (^ {i) J 

then the process X^^^ attains the point b before a.s., for any x £ (0,6). In that case, the 

process never hits a.s., because continuity in time requires that it pass through the interval 

(0, b) to reach and the Markov property requires that each time it re-enters the interval it 

must exit through b. The same statement holds for the right endpoint: if 

then the process X^^r attains the point b before reaching oo a.s., for any x £ (6, oo), and thus 
never reaches oo in any finite time a.s. We note in general that 



[-l^^f]=C[a\x)P{x)]-\ 



il>{x; b) := exp{- -^dt } = C [a'^{x)P{x)\ , (4.38) 



using ( [4.35| ), where the constant C = a'^{b)P{b). Thus, for any of the processes ( [4.24 ), ( 4.25| ), 
it follows that 

^(x;b) = C'—^^ (4.39) 
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using o"^(x) = 2x and the formula ( [4.29| ) for the invariant measure. Clearly, L^ = L xp{x; b) dx 
and L^ = f^ Tp{x; b) dx. 

In the case at hand, for the equipartition solution, we see that ipQ{x;b) = Cx~^e^. Since 
c > 1, both L^ = +00 and Lf = +cx). Thus, we conclude that the process never reaches the 
boundary in finite time a.s. This completes the convergence proof for the case of bounded A{0). 
By itself, this is only a strengthening of the result in section 4.1, since ^(0) € L°° implies that 



A{0) is L^ with respect to the weight wdx) 



-,c— 1 _— x 



However, the sense of convergence is 



stronger, being now pointwise in x uniformly on compacts. 

Having proved the result for bounded ^(0), we now extend to the case where r(0) G L'^ 
and r(0) € L^ with respect to weight x'^~^. We have, using r(a;,0) = yl(x,0)e~^, that 



|A(x,0)|<||r(0)||Lcce^ 



(4.40) 



and 



00 1 roa 

\A{x,0)\Po{x) dx - 
J^ ic 



1 f°° 

-- \T(x,0)\ x^-^dx < 00. (4.41) 

(c) Jo 

Because of the latter result, we may choose M > 1 so large that for any small e > 0, 



\A{x,0)\Po{x) dx< e. 



(4.42) 



M 



Let us then take 



A^^\x,0):-- 



A{x,0) x<M 
x> M 

Thus, by (glO| ), A^^^\0) is bounded: ||A(*^)(0)||l- < ||r(0)||L-e^^. By triangle inequality 



(4.43) 



dy Aiy,0)P{y,T\x,0) 



dy A{y,0)Po{y) 



< 



dy A'^^'\y,0)P{y,r\x,0) - / dy A'^''\y,0)Po{y) 



+ 



M 



dy |>l(y,0)|P(y,r|x,0) + / dy |^(y, 0)|Po(y). (4.44) 

JM 



To control the second term we employ the estimate 

P(y,r|x,0) 



sup 



< 



y^o Poiy) -{l-e-r 



(4.45) 
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This is proved using the inequality ([p|], 7.3.2(4) and 7.2.2(12)): 



„ I Re 2 

\z-''Uz)\ < (4.46) 



r\l/2 



with z = ^ ,_ _/ — and u = c— 1 in ( [4.33| ) . Thus 



P{y,T\x,Q)< ^^_^_,^, exp| -^-^ 1. (4.47) 

The maximum of the exponent is found to occur at y = xe^ , yielding the estimate ( [4.45| ). That 
inequality can then be used to bound 

r dy \A{y, 0)|P(y, t\x, 0) < ''"_ e (4.48) 

Jm (1 -e ^r 



because of the condition ( 4.42| ) on M. Since ^(*^^(0) is bounded for a fixed M, our preliminary 



convergence result applies and we obtain 



lim sup 



< (e^ + l)e. (4.49) 



/ dy A{y, 0)Piy, t\x, 0) - / dy A{y, O)Po(y) 
/o Jo 

Since e was arbitrary, we get 

lim/ dy A{y,0)Piy,T\x,0) = dy A{y,0)Po{y) := cq. (4.50) 

■^^°° Jo Jo 

and therefore 

lim r(a;,r) = co$o(a:) (4.51) 

r— >oo 

uniformly on compact subsets of (0, cxd). 

To identify the constant cq , we note that the Corrsin invariant of the initial data is 

K{0)=Y~^uJd-i A{x,0)<^o{x) x""-^ dx. (4.52) 

Furthermore, Po{x) = ■^j'^^^LOd-ix^~^^o{x), where K is the Corrsin invariant of the equipar- 



tition solution <l>n- Thus, 



K{0) = K ■ A{x, 0)Po{x) dx = K ■ cq. (4.53) 

Jo 
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Of course, this is consistent with the tinie-invariance of the Corrsin integral, which, by (4.51), 
gives K{t) = cqK for all r > 0. 

If the Corrsin invariant is finite but vanishes, then the results in the preceding Section 4.2 
suggest that the asymptotic behavior will be described by ^p, if Jp is the first non- vanishing 
invariant. We next prove such a result, which is relevant to initial conditions with power-law 
decay for z^ > 7. The theorem we establish is the following: Suppose that 0(0) G L°° and 
J^ddr rP'^\@{r,0)\ < 00. If Ji = 0, i = 0, l,...,p — 1 but Jp 7^ 0, then the solution rescaled 
appropriate to parameter a = d + p'j satisfies 

lim r(„)(x,r) = Cp%{x) (4.54) 

uniformly on compact subsets of x. Furthermore, the constant Cp is the ratio Jp/Jp, where Jp 
is the value of the invariant for the scaling solution Qp(r,t) = 'd'^{t)^p{r/L{t)). 

We give the proof first for the bounded A(0). We wish to make a proof very similar to the 
previous one. In fact, 

A{x,t) = E[A{X,^^M (4.55) 

where X^^t is the diffusion process appropriate to a = d + p^- It is enough to show that 

/■oo 

lim A{x,t) = / A{x,Q)Pp{x) dx := Cp (4.56) 

uniformly on compact sets of x with 

Ppi^) = f^^^ [L^-Hxfe-^x'^-' (4.57) 

the invariant measure of the process. The complication is that the diffusion in all the cases 
p > 1 has singular points and decomposes into p+ 1 simple pieces, each with its own invariant 
measure. In fact, using ( [4.24 ) for the drift a{x) and ( ^.261) for <I>p(x), 



" 2x 



a{x)=c-x + Y,^^— (4.58) 



fc=i ^ ^'^ 



where x^, k = 1, ...,p are the p roots of the generalized Laguerre polynomial L^ ^{x). We recall 
that these are all real and simple, and located in the interior of the interval (0,oo). Cf. |49|, 
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Theorem 3.3.1. We may label them in increasing order xi < X2 < ■ ■ ■ < Xp. Because of the pole 
terms, the zeros are repulsive singularities of the drift field a(x) and the regular set of points 
of the process decomposes into a disjoint union of p + 1 open intervals, Ik = ixk,Xk+i), k = 
0, 1, ..,p, with xo = and Xp+i = oo. It is easy to calculate that for each x,b € I^, 

with C some constant. Thus, /^ il){x;b) dx = +oo and /^ ^^^ ip{x;b) dx = +oo, so that, again 
by [O], Section 21, Theorem 1, the boundary points of each interval are inaccessible in finite 
time a.s. Thus, the process is not ergodic but instead there exist p+\ distinct, ergodic invariant 
distributions Pp supported on the intervals /fc, k = 0,1,..., p. Up to a normalization factor 
Wk = {Jj Pp{x)dx)^^ , these coincide with Pp restricted to the interval Ik, i.e. 

Pj,'\x) = WkPpix)\,^ (4.60) 

On each interval separately, the Theorem 3, Section 23 of |5l| applies. Thus, for x £ Ik 

lim E [A{X^^r,0)] = Wk f A{x, 0)Pp{x) := Cp^k- (4.61) 

We need to show that the constants Cp^k are, in fact, independent of k. 

By assumption, the initial data has the first p invariants vanishing, Jq = ■ ■ ■ = Jp-i = 

but Jp / 0: 

^/' , Aix,0)L'f-\x)L^-\x)e-''x^-Ux = (4.62) 

^{P + c)Jo ^ 

for ^ = 0, ...,p — 1, and 

p\ 



/■oo 

The uniform bound |y4(x,r)l < ||y4(0)||ioo follows from ( [4.55| ). Thus, the J^, £ = 0,1,. ..p are 
rigorously dynamical invariants and the equations (4.62|),(p.63D hold with A(x,0) replaced by 



A[x,t). Then we may apply ( [4.61| ) and Lebesgue's theorem to conclude 

^^— Y^ Cp^k I '^' Lf^{x)Ll-^{x)e-''x^-^ dx = (4.64) 



r(p 



A;=0 ^'^ 
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p\ J-, r^k+1 



for £ = 0, ...,p — 1, and 

Tini a) T.^P''^ I ^^^ [LrHx)]'e-^x'^-' dx = J,. (4.65) 

These may be summarized as a matrix equation Mcp = J with 

M^,fc:= / L''f^{x)Ll-^{x)e-''x''-^ dx. (4.66) 

By orthogonahty of Laguerre polynomials one solution is Cp^k = Jp/ Jp for all k, where recall 

1 f^ 

Jp := :p7^ / ^p{x)L'l'^{x)x''~^dx 
r(c) ^0 ^ 

is the pth invariant of the scaling solution ^p itself. This is the unique solution if the matrix M 
is nonsingular. That is equivalent to the statement that an arbitrary polynomial p{x) of degree 
p can satisfy 

/ p{x)Ll'^{x)x''-^dx = Q, k = 0,l,...,p (4.68) 

only if p(x) = 0. In fact, this is true, because the polynomial has at most p real roots but there 
are p + 1 intervals. Hence, there is at least one interval Ik, k = 0,1, ...,p on which it does not 



change sign. Then (4.68) implies that p{x) = on that interval, and, hence, everywhere. Thus, 
we conclude. As a by-product of this argument, we have shown that 

/•OO -1 /"OO 

Cp= A{x,0)Pp{x) dx = --- T{x,0)L'p'^{x) x'^-^dx (4.69) 

Jo r(cj Jo 

is given by Cp = Jp/ Jp, as claimed. 

Having proved the result for bounded ^(0), we now extend to the case where r(0) E L°° 
and r(0) E L^ with respect to weight x'p~^^~'^. Here we make use of the observation that the 
scalings by ap = d + p7 and ao = d are simply related by r(„ ){x,t) = eP'^r(^)(x, r), where we 
have made the a-dependence explicit. After this we shall employ the equipartition scaling and 
the corresponding definition of r(^)(x,r) = j4(^)(x, t)e~^, so, when no a-dependence is given, 
a = d is implied. Then, we must show that 

lim eP^A{x,T) = Cp^Ll-\x) (4.70) 



T^OO 
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uniformly on compacts, where 

1 r°° 

cp = ,^^ T{x,^)Ll-\x) x'^-^dx. (4.71) 

We define the following auxilliary function: 

/•OO /'OO /'OO 

H{v, 0) := (-l)P / dyi dy2--- dy^ y^-^r(yp, 0). (4.72) 

Jy Jyi Jyp-i 

It is not hard to see that reversing orders of integrations by the Tonelli theorem gives 

/ \H{y,0)\dy < / dyp dyp-i • • • / dy y;-^\Tiyp,0)\ 
Jo Jo Jo Jo 

= - / \T{y,0)\yP+'-'dy<oc. (4.73) 

p! Jo 

Of course, it follows directly from the definition that 

1 dP 

Because of the vanishing of the first p invariants, one may readily check that 

d'H 



dy' 



-(y,o) 



= 0, £ = 0,l,...,p-l. (4.75) 

j/=0 



(4.76) 



Thus, making p integrations by parts with -^L^ ^{y) = ( — 1)^ 

— -/ H{y,^)dy = —- H{y,0)L;-\y) dy 

r(c) Jo r(c) Jo dyP 

^ ml r(y,0)L^-i(y)y^-idy = Cp. 
Since the integral is absolutely convergent, we may choose M sufficiently large that 

-1^ r\Hiy,0)\dy<e. (4.77) 

r(c) Jm 

Let us define a decomposition H{y, 0) = i?^*^) (y, 0)+;^ ^^^^ (y, 0) via i?^*^) (y, 0) = c/?^*^) (y)i?(y, 0) 
and H {y, 0) = Tp^ >{y)H{y, 0) for a smooth decomposition of unity tp^ '{y) +Tp^ \y) = 1, 
with (/7(*^)(y),^(*^)(y) > and <p^'^^\y) = for y > M + 1, ^^^'^\y) = for < y < M. We 
may then define a corresponding decomposition r(y,0) = r(*^)(y,0) + T (y,0) via ( 4.74| ), 
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and likewise for A{y,0) = A^^'^>{y,0) + A (y,0) and Cp = Cp + Cp . Next we can employ 
the transition probability of the a = d process in a triangle inequality: 



pPT 



/ dy A{y,0)P{y,T\x,0) - —^L;-\x) dy H{y,0) 

Jo l[C-\-p) Jo 



< 



r{c + p) 

zP^ / dy AW(y,0)P(y,r|x,0) - cW-^l;-1(x) 



/o 



r dy /'''\y,0)Piy,T\x,0) 
Jm 



+ f(^l^^ '(^) 



dy 



M 



H'''\y,o) 



(4.78) 



As before, p(*^)(0)llL°° < ||r(0)||L-'e(*^+^). Thus, we can identify the first term as 



Ai':\x,r) 



.im 



■7pp|Lp~ (x)| and appeal to our preliminary result for L°° initial data to conclude that this goes 
to zero uniformly on compact sets of x. Of course, the third term is less than -7^\Lp'^^{x)\ ■ e 
by the assumption ( 4.77| ). The main problem is to control the middle term. 

If we substitute the expression ( [4.74| ) for T [y, 0) in terms H (y, 0) and integrate by 
parts p times, we obtain 

dV\p{y,T\x,oy 



dy A^^'\y,0)P{y,T\x,0) 



dyH^^'\y,0) 



iM JM \ dy 

Employing the formula ( [4.33 ) for the transition probability gives 

rdyA^''\y,0)P{y,T\x,0) = -'^''~' 
Jm 



yC-lg-J/ 



(4.79) 



(1-e-^)^ 



M 



..F-(.o) l-^J 



exp 



X + y 



Z-(^-l)/e-l(^) 

(4.80) 



with z = ^ ,_ _l — . The derivative can be evaluated by the generalized product rule DP{uv) 



X]r-n ^ r) D'^u ■ D'P '^v and the relation 



jr=o yrj 



2xe 



d 



We note that 



_(f_ 



■exp 



d 



dy (1 — e '^y zdz 



■^ ' y { I -i\r I r -i\—r I X + y 

(—1) (e — 1) ■ exp 



Likewise, defining ^ = x/{l — e 
d- 



dy^ 



z-(^-i)4_i(z) 






' z'^'+^-^h,^r-l{z) 



(4.81) 



(4.82) 



(4.83) 
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using (J^Ii\ ) and {^Y [z-^I^iz]] = z-(^+'') /;,+,. (z), |2|, 7.11(20). Summing over all the con- 
tributions and using the estimate ( [4.46| ) for the Bessel function gives 



dP 
dyP 



expi-^U-(-^)/e-i(.) 



< 



■ exp < 



1/2' 



2c-ir(c + p)(e^ - 1)P 



1 



We employed the relation 



L'-Hx) 



E 

r=0 



c+p-1 
p—r 



i-xY 



ip] 



(4.84) 



(4.85) 



in |25|, 10.12(7). Thus, |Lp(x)| < Lp ^(—x). The exponential term is the same as was shown 
before to be bounded over all y by e^. Thus, we obtain finally from (4.80) that 



pr 



e- I dy A^^\y,0)P{y,T\x,0) 

M 



< 



[i-e-^Y+p r{c + p) 



M 



dy\H{y,0)\ (4.^ 



Thus, 



lim sup 



-pr 



AM) 



dyA'"''{y,Q)P{y,T\x,Q) 



M 



{c)p " ^ ' ' 



(4.87) 



by the definition ( 4.77 ) of M. Since ^(*^'(0) is bounded for a fixed M, our preliminary conver- 
gence result applies and we obtain 



eP^ I dy Aiy,0)P{y,r\x,0)-Cp^L;-\x] 

ICjp 



lim sup 

T— >00 

Since e is arbitrary, we get 



<(e- + l)^L-i(-x).6. (4.^ 
(c)p 



lim eP'^A(x,T) = Op-f—Vl ^(x) 

T^OC ''(C)p P 

uniformly on compact subsets of (0, oo) and likewise 



(4.89) 



lim Tap{x,T) = Cp^p{x). 



(4.90) 



The most surprising consequence, in view of traditional beliefs, is the implied breakdown of 
spatial PLE on length-scales comparable to the integral scale L(t). Naively one would expect 
that PLE holds in space for all long-range powers laws, ^ Ar~" with any a > 0. Instead, 
the memory of the initial conditions is through the invariant Jp and not through the initial 
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amplitude A of the power-law tail, when a > d + p^ and Jp ^ 0. Of course, one expects that 
PLE still holds on the logarithmically larger length-scale L^,{t). We shall consider that below. 
We caution that the above are pure PDE results, and do not necessarily describe statistically 
realizable situations. As already discussed in Section 4.2, it is not possible that Jq = Ji = 
but Jp 7^ 0, oo for some p > 1 is consistent with realizability. Thus, only the p = 0, 1 cases of 
the above theorems are relevant to solutions arising from the Kraichnan model. If realizable 
initial data have a long-range power q > d + 7 and if Jq = 0, then it must be that Ji ^ 0, so 
the second theorem with p = 1 applies and the appropriately scaled solution converges to $1. 
This agrees with the physical picture presented in Section 3. 

It would be interesting to prove spectral analogues of these theorems. We shall not attempt 
to do so here. However, we note that the conditions of the present theorems can be implied 
by spectral conditions on initial data. In fact, an initial condition with power-law spectrum 
~ Ak°'~^ for an exponent a = d+u, i^ > 7 has spatial decay r~° at large r (except when u = 2j 
an even integer, when the decay may be even faster). Thus, Jp ^ 00 for p = 0,1. However, 
Jo = 0, because the Corrsin invariant is proportional to the coefficient of the k'^~^ term of the 
low-wavenumber spectrum, which is assumed to vanish. Thus, the conditions of the second 
theorem apply, for Ji ^ 0. In that case, the long-time limit for the scaled scalar spectrum 
F{K,T{t)) := E{K/L{t),t)/{^'^{t)L{t)) is presumably ciFi(k) with ci = Ji/Ji. Although the 
initial spectrum E{k, tg) had the low-wavenumber form ~ Ak°^~^ for k somewhat smaller than 
-^^^(^o)) at sufficiently long times the spectrum is ~ cA{i',d)L''[t)k'^~^^ for wavenumbers k 
somewhat smaller than L~^{t), with the same constant A{v,d) as in (|3.68D . 
(4-4) Initial Data with Slow Decay & No Finite Invariant 

The last case to consider is that all of the invariants Jg, I = 0,1,2,... of the initial data 
r(x,0) are either zero or infinite. It is important to appreciate that this is the case for the 
nonexceptional scaling solutions ^ themselves. In fact, for pj < v < {p + 1)7, formula (3.40) 



implies that ^ has the invariants J^ = 0, £ = 0, ...,p — 1 and Jg = ±00 for all i > p. Because 
these are conserved by the dynamics, only initial data with the same pattern of invariants can 
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exhibit dynamical scaling at large times, governed by a nonexceptional ^. From this point of 
view the failure of PLE that was discussed in the last section is not so surprising. To observe 
a long-range power-law in the scaling limit, the initial data must have the same invariants as 
the final scaling solution. 

We prove now some theorems which establish convergence to the nonexceptional scaling 
solutions. We consider first the easiest case, when —d < v < 0. We show the following: 
Suppose that 0(0) G L°° and that with < a < d 

G(r, 0) ^ Ar-'^, r » Lq. (4.91) 

Then, i/r(Q)(x,r) is the solution rescaled corresponding to the parameter a, 

lim ?(„) (x, r) = Cc,$(„) (x) (4.92) 

T— >00 ^ ' ^ ' 

uniformly on compact sets. The constant Ca = A/ A where A = 7" — t^jy is the constant 



prefactor in the asymptotic power-law ( 3. 23( ) for the scaling solution ^ia){p) = ^ ( ~) i ~~ ) • 



Thus, spatial PLE holds. This seems to be the optimal result that could be expected. The 
proof depends upon the observation that r(x,0) may be written in the form 

r(x,0)=yl(x,0)<^(„)(x) (4.93) 

where ^(0) € L°° and lima;_>oo A{x, 0) = c^ < cxo. In fact, it follows from the integral represen- 
tation (|3.45| ) that <I>(q)(x) has no zeros and decreases monotonically from the value 1 at x = 0. 
This, along with asymptotic power-law behavior which matches that in the initial data, gives 
the existence of a multiplicative perturbation A{x, 0) with the required properties. 



The proof uses the stochastic representation by the diffusion with drift a(x) given by ( 4.25 ) 
and diffusion 6(x) by ( g^ej) . By (|3l5|) 



^'(a) (^) fn ue-^^n'^-i (1 - uY-^-^du 
^ > = -do !^ i g („i 0) (4 94) 
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Furthermore, ^'"^ >, ~ — ^ for x » 1 by ( 3.22| ). Thus, we see that \a{x) — c\ < x and that for 



x> 1 

a{x) r^c-2a + x + 0{x-^) (4.95) 

Therefore, the drift is nonsingular except at the right boundary x = oo of the interval, which 
is infinitely attractive. The left boundary x = is repulsive. For any x, b in (0, oo) we 



have 'ilj{x;b) = C "i'i ", . by ( [4.39| ). Because ip{x;b) r^ C'x^" for x ^ 0+ and c > 1, L, 



/q ip{x;b)dx = +00 and the process Xx^t never reaches a.s. However, il){x;b) ^ C"x'^°'~^e~^ 
as X ^ +00, and thus L^ = J^ '4){x;b) dx < +oo. In that case, [^], Section 16, Theorem 1 
implies that 

lim Xx^T = +00 a.s. (4.96) 



T— >00 



Since ^4(0) E L°° and lim^^oo ^{x, 0) = c^, we may apply Lebesgue's theorem to conclude that 

lim A(x,r) = lim E [^(X^,^,0)] = c^. (4.97) 

r — ^oo T — >oo 

The proof is complete. 

We now consider the case v > Q, v ^ ^7, £ = 0, 1, 2, .... This case is more difficult because 
the Kummer functions <l>(a, c;— x) with a > c have positive, real zeros. This is, of course, 
required by the vanishing of a certain number of the J-invariants. However, it presents a 
difficulty for the analysis. We prove the following: Suppose that the initial data B(r, 0) is a 
bounded perturbation of ^(a) /o*" a > d, a y^ d + £7, ^ = 0, 1, 2, .. in the sense that 

r(p,0)=^(p,0)<&(,)(p), (198) 

with A{0) S L°° and that limp_+oo A{p, 0) = Cp < 00. Suppose also that Jq = Ji = • • • = Jp-i = 
0. Then, 

lim r(c,)(x,r) = Cp$(a)(x) (4.99) 

uniformly on compact sets. The conditions on A{0) imply that r(p, 0) ~ Ap~°' for p ^ Lq, 
with A = CpA as before. Thus, spatial PLE holds for this class of initial data. However, the 
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result is not optimal, because the conditions also imply that r(p, 0) has zeros at precisely the 
locations of the zeros of ^(a){p)- The condition on vanishing invariants — which is necessary for 
the convergence result to be valid — implies the existence of such zeros, but not at precisely the 
location of the zeros of ^(a)- There will be more general initial data in the domain of attraction 
of the scaling solution, but there is then a difficult problem in accounting for the motion of the 
zeros in time. We leave that to future work. 

The proof of the above result uses again the stochastic representation. We note that the 
Kummer function ^{a, c; —z) is an entire function in the complex z-plane and, for real a, c has a 
finite number Z of zeros, all simple (|2^, Section 16). It then follows from Weierstrass' product 
formula for entire functions that ^{a, c; —z) = e^'^' Y[k=i(^~'^k) where g{z) is an entire function 
and the product is over the Z complex zeros. When a, c > the real zeros of ^{a, c; —x) are all 
positive and there are precisely p such positive, real zeros xi, ...,Xp when c-\-p<a<c + p + l 
(11, Section 16). Thus, 



$; Jx) P 1 

where h{x) is a C°° function on the real line. For x ^ 1, ^ (^\ ~ — f because of the asymptotic 
power-law form of ^la)- Thus, we see that the drift field is 

P 2x 
a(x) = c + x + y^ \-2xh{x) (4.101) 

^-^ X — Xk 

fc=l " 

with p repulsive singular points at the zeros of <!>(„) (x) and, as before, 

a(x) ~c-2a + x + 0(2;"^) (4.102) 

for X ^ 1. Thus, the process decomposes into p + 1 simple pieces and the regular set of 
points consists of the disjoint intervals Ik = (x^jX^-i-i), with xq = 0,Xp = c«. For each of the 
first p intervals, L^ = Lf = +oo as before. The process is ergodic on these intervals with 
invariant measure Pk{x) = WkP[a){x) for P(a){x) = 3;'^~^e^'5L)(2;) and w^ = fj P(q)(x) dx 
However, on the final interval Lj~ = +(X) but Lf < +oo. Thus, liniT-^oo ^x,t = +oo a.s. when 
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X (z Ip. Therefore, we have that 

hm E [A{X^^r, 0)] =Wk [ A{y, 0)P(«) (y) dy := Ck (4.103) 

when X (z Ik, k = 0, ...,p — 1 but that 

hm E [A{X^,r, 0)] = A{+oo, 0) = Cp (4.104) 

T—tOO 

when X £ Ip. We must show that c^ = Cp for A; = 0,l,...,p — 1. Using the above convergence 

result and the vanishing of the first p invariants for the initial data, one can show as before that 

the c-coefficients satisfy 

p ,. 

^Ck I Ll-^{x)^(,^){x)x^-'^dx = Q (4.105) 

for ^ = 0, ...,p — 1. This can be written as an equation for the vector c = (cq, ...,Cp_i)^ of the 
form Mc = d with 

Mtk= I Ll-^{x)^(^a){x)x^''^dx (4.106) 

for i,k = 0, ...,p — 1 and 

de = -Cp f L^"^(z)$(„)(x)a;"-^ dx (4.107) 



for i = 0,...,p — 1. Because the first p invariants vanish also for <1>(-q,)(x), one solution is 
Co = • • • = Cp-i = Cp. In fact, the matrix M is non-singular, because the same argument as 
before shows that a polynomial p{x) of degree p — 1 which satisfies 



p(a;)$(„) {x)x''-' dx = (4.108) 

for k = 0, ...,p — 1 must vanish identically, p{x) = 0. Thus, the solution is unique, as required. 
The previous results again support the physical picture proposed in Section 3. As before, 
the theorems of this section do not describe a statistically realizable situation when a > d + 7, 
although, as pure PDE results, they are valid. When the initial data has a long-range power 
decay with exponent a < d + j, then the theorems for p = 0, 1 do apply. In that case, two 
possibilities occur depending upon the value of Jq. The additional conserved quantity Ji plays 
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no role because, for a < d + 7, Ji = cxo always. If < Jq < 00, then the long-time behavior is 
that associated to ^q, as shown in the previous section. However, if Jq = or 00 (the latter will 
always hold for a < d), then the theorems of this section imply that the long-time behavior is 
that associated to the realizable scaling solution ^ta) for the same a. This is the circumstance 
in which spatial PLE holds, in the conventional sense. 
(4-5) View on a Larger Length-Scale 

Finally, we will discuss briefly the scaling by L*(t) logarithmically bigger than L{t). Our 
analysis will be only heuristic and no proofs given. The solution to the equation ( |3.6D is 
L{t) = jDi{t — tg)^ foi' a virtual time origin t^ related to the initial data Lq. Let us define 



L,it) = jDi[{t - t*o)log{t - t*o)] 



*\]i/'r 



(4.109) 



We now study rescaled solutions 



e{r,t) = ^\t)T,{r/L,{t),T{t)), 



(4.110) 



where i?(t), T(t) are as before. In particular, T{t) := log(t — tg) = Lj{t)/L'^{t). The function F* 



solves: 



' dT 






1 1 d 



dp 



P'^'-'^-§fiP^-) 



where p = r/L^{t). If we make the change of variables x = 7 p'^, we obtain 



+ (l + ^)p^(p,r)+ar,(/,,r). 



^-— (x,r) =x——{x,t) + aT:,{x,T) +- 

OT OX T 



X „ * (x, t) + {c + x)-— -^(x, r) 



dx'^ 



dx 



(4.111) 



(4.112) 



with a = a/7, c = d/7. 

If the initial datum has the power-law form r(x, 0) ~ Ax~"' at large x, then we expect that 
r^,(x,r) converges to that same power pointwise in x. This motivates us to define A^{x,t) via 



r*(x,r) = X "'A^{x,t). 



(4.113) 



Thus, we want to show that liraT^oo A^{x,t) = A[-\-oo,0). Substitution of ( 4.113 ) into ( [4.112| ) 
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yields the equation 



dA dV 1 



9^r* 5r* aia — c + l — x) 

X „ „ [X, r) + (c - 2a + x)——[x, t) -\ A^{x, t] 



dx"^ 



dx 



X 



(4.114) 
Heuristically, we may expect that, at long times, the terms in the bracket can be neglected 
because of the factor 1/r. In that case, the solution is governed by just the hyperbolic equation 
-g:^(x, r) = x-Q^{x, r). This is solved by the method of characteristics, as A^{x, r) = A{e'^x, 0). 
If this crude approximation is valid, then liniT-^oo ^♦(a^)''") = ^(+oo,0), as required. 

A more refined estimate can be made by using a stochastic representation corresponding to 
the diffusion process with time-dependent drift 

c — 2a + X 



a(x,T) = X + 



(4.115) 



and diffusion 



b{x, t) = cr2(x, r) = — . (4.116) 

T 

If Xx^T is the process started at x at time tq, the solution of ( [4.114 ) has the representation 



A^{x,t) 



E 



exp ( / c{Xx^„,a)da] A{Xx^r,0) 

To 



where 



c{x,t) 



a(a — c+ 1 — x) 



X ■ T 



(4.117) 



(4.118) 



See [^], Section 11, Theorem 3. Because the drift a{x,T) ~ x at large x,t and, furthermore, 
o-^(x)/x^ — > as X ^ oo or r — > oo, one expects that X^^r ~ e"^ as r — > oo a.s. See |51|, Section 
17, Corollary 2. The contribution to the integral of c(x, cj) in the exponent for a^r may be 
estimated by c{X,j.^a,(^) ^ —a/a when X^^^ is large, as it is with high probability when a^T. 
This gives a term ~ exp(— alogr) = r^", which vanishes for large r. On the other hand, the 
contribution from ct^tq is of the order of exp (0(1)/t) ~ 1 for large r. Thus, we expect again 
that limr-^oo A:^{x,t) = ^(+oo,0). 
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5 Conclusions 

In this paper we have studied in detail the 2-point correlation function in the Kraichnan model, 
which, for homogeneous and isotropic conditions, is governed by the equation (|1.9D. Our main 
technical results are as follows: 

(i) We have found all the possible self-similar decay solutions, parameterized by space di- 
mension d, velocity field roughness exponent (^ = 2 — 7, and a scaling exponent a > 0, in 
terms of Kummer confluent hyper geometric functions. 

(ii) We have shown that these solutions are statistically realizable for all a in the interval 
(0, d + j] and for no a in the complement (d + 7, 00). 

(iii) We have exhibited an infinite sequence of invariants Ji, i = 0,1,2,... of the equation 
( |l.g| ), defined by the formula ( [4.15| ). Jq coincides with the well-known Corrsin invariant. 
At least one of Jq, Ji must be nonzero (possibly infinite) for realizable initial data. 

(iv) In terms of these invariants, we have identified initial data in domains of attraction of the 
scaling solutions found in (i). Our theorems relevant to the Kraichnan model covered the 
following cases: 

Jg = Ji = 00 (r(a)(T") -^ ^{a), < a < d) 

Jo = 0, Ji = 00 (r(Q^) (r) -^ $(„) , d< a < d + ^) 

< Jo < 00, Ji arbitrary (r(Q)(r) -^ ^q, a = d) 

Jo = 0,0 < Ji < cx) (r(Q)(r) ^ $1, a = d-h7) 

The first two results apply for initial data with a power-law decay ^ r~°' for a < d + 'j, 
whereas the last two include (among other possibilities) initial power-laws with a > d + 'j. 

There are several questions left open in this work that deserve to be addressed. It would be 
worthwhile to establish convergence for more general initial data with slow decay and no finite 
invariants, allowing for a non-coincidence of its zeros with those of the final scaling solution. 
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It would be interesting as well to extend the convergence results to the spectral domain and 
to establish convergence to a power-law solution on the logarithmically larger scale. Perhaps 
the most interesting issues not treated in the present work are to consider scalar decay in 
the Kraichnan model with initial anisotropy and/or inhomogeneity and to consider the self- 
similarity (or not) of the higher-order A^-point scalar correlations with N > 2. The latter is 
relevant for realistic decay problems, where no closure of the moment hierarchy is found. In 
fact, von Karman and Howarth in their classic analysis [^ made a self-similarity hypothesis for 
the 3rd-order as well as the 2nd-order velocity correlations, in order to derive a scaling equation. 
Perhaps the most novel result in this work is the discovery of the infinite sequence of invari- 
ants Ji, £ = 0, 1, 2, ... in the Kraichnan model. It would be most interesting to know whether 
there are analogues of such invariants for true passive scalars, besides the Corrsin invariant 
K = Jo, or, for that matter, for Navier-Stokes turbulence. In the latter case, the integral in- 
variant C = Jjgd c?^r BllI^t) found by Saffman [£2[ is analogous to the Corrsin invariant, while 
the integral A = J^^ d'^r r'^B^L^r) found by Loitsyanskii [£^] is rather analogous to our next in- 
variant Ji. However, it is well-known that A(t) is actually time-dependent, due to the effects of 



spatial long-range correlations [29|. The connection of the "generalized invariants" introduced 
here with the "zero modes" studied in ^^ is very intriguing, especially as the latter exist also 
for anisotropic statistics and, plausibly, for higher-order A^-point correlators. 

Our results confirm in the Kraichnan model the "two-scale" decay picture earlier found for 
Burgers equation in p|. The fact that there are no realizable scaling solutions for a > d + j, 
implies that initial data with such a r~° decay at large distance cannot exhibit a self-similar 
decay with the same power-law behavior on the scale of the integral length L{t). Instead, we 
have shown that a wide class of such initial data converge on that length-scale to the self- 
similar decay solution Qi{r, t) corresponding to a = d + 7. A heuristic argument supports the 
conclusion that the initial power-law decay will be preserved on a length-scale L^{t), which is 
logarithmically bigger than L{t). The fact that the same phenomenon occurs for two such very 
different models as Burgers equation and Kraichnan's passive scalar model argues for its gener- 
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ality. In fact, the physics behind it appears to be just the phenomenon of "backtransfer" , which 
closure results indicate is a very general feature of turbulent decay. The specific, simplifying 
features of the models only enter in calculating the rates of this backtransfer. 

We therefore expect that the results of this work will have relevance for real passive scalars 
and other turbulent decay problems. It is thus of some interest to specialize our results to 
that case of the Kraichnan model which most closely mimicks the true passive scalar, namely, 
d = 3 and (" = g. In terms of the scalar low-wavenumber spectral exponent n = a — 1, there 
is a self-similar decay with PLE if initially — 1 < n < |. However, if | < n < 4 then these 
traditional expectations are violated. In the whole range "- > | the scalar energy decay law 
is E{t) ~ (t — to)"^"'^ J as traditionally expected only for the lowest value n = |. We have 
learned of unpublished work of D. Thomson who finds exactly the same behavior in a simple 



model of a mandoline source [47, ^]. This leads us to believe that the results presented here 



have some validity outside the white-noise model. In fact, on this basis, we have presented a 



phenomenological extension of this picture to decay of 3D incompressible fluid turbulence 1 54 1 . 

Although the physics will be different, similar phenomena may occur in other nonequilibrium 

decay processes. During the preparation of this paper we became aware of a "two-scale" picture 

for phase-ordering by Cahn-Hilliard dynamics with power-law correlated (or fractally clustered) 



initial data |55|. From a renormalization group (RG) point of view, our results correspond to 
a case where some of the "fixed points", i.e. the scaling solutions, are unphysical. A careful 
study of this example using general methodologies such as RG may help to illuminate subtleties 
in their application. 
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APPENDIX: 

SELF-SIMILAR SCALAR SPECTRA FOR A BROWNIAN VELOCITY FIELD 

In general, it is not possible to evaluate explicitly the Fourier transforms of the scaling func- 
tions, ^{k), not even in terms of standard special functions. This is not possible even for the 
equipartition case q = d, in which case the Fourier transform is the spherically symmetric Levy 
stable distribution with parameter 7 in d dimensions. However, it is known that the stable 
distributions are calculable explicitly for the parameter value 7 = 1, in which case they are 
the d- dimensional Cauchy distributions. Since then also C = 1, this case corresponds to an 
advecting velocity field with Holder exponent H = ^, 01 a Brownian-type random field. 

In fact, the Cauchy distribution of parameter f3 in d-dimensions is defined by the Fourier 
transform 

with the result 

^^"'^) = (i^o^F(iy • (^2+,2)i(.+i)- (^•2) 

This follows from |25|, 7. 7.3(16) with the parameter choices 7 = /?, a = k, /x = i(d — 2),p = 



2(^ + 2) there, because the hypergeometric series then simplifies to a binomial series for the 
function ( 1 + §7 ) • We shall exploit this fact to evaluate the Fourier transforms ^{k) 

for 7 = 1, at least for all of the exceptional cases a = d + i, i = 0,1,2, .... We may use G{k; (3) 
as a generating function in those cases, because 






(5.3) 



/3=1 



as a consequence of ( p. 26 ) 



The results obtained by use of this formula are, explicitly, for ^ = 0, 



^ i -. 2 r(d) 1 , ^, 
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for 



for 



for 



1, 



2, 



3, 



cDi(k) 



^2{k) 



T{d) 



(^+1)^2 



^3{k) 



2 r(d) (d + 3)[(d + 2)^6 -3^4] 



(47r)'^/2 r (f ) ci(d + 2)(l + K2)^('^+7) ' 



(5.5) 



(5.6) 



(5.7) 



and for 



$4(k;) 



2 r(d) (d + 4)((i + 2)K8-6(d + 4)K6^3K^ 



(4^^)d/2Yfd^ 



,2^|{rf+9) 



(5.8) 



,2/ (i((i + 2)(l + K2 

Of course, these transforms and those for ah i are analytic at k = with radius of convergence 
of the power series equal to 1. It may be easily checked that the coefficients are equal to the 



Bj given by (3.37). In particular, there is no singular contribution at k = 0. 

These explicit results illustrate several of the general conclusions reached in the text. For 
example, we see that ^o{k),^i{k) are everywhere nonnegative and that <I>i(k) = ~]i'^^;^{''^)- 
Furthermore, in agreement with the result that realizability fails whenever i/ > 7, we see that 
the Fourier transforms for i = 2,3,4 are not everywhere positive. In the cases ^ = 2,3 the 
coefficient of the dominant low-wavenumber power is negative, so realizability fails for the 
lowest wavenumber range. However, for i = 4, realizability fails despite the coefficient of the 
dominant low-wavenumber power being positive. In fact, the polynomial 



P(k^) = {d + 4){d + 2)k^ - 6{d + 4)^6 + 3k^ 



(5.9) 



has two positive roots k^ = -^^ 



1 + / 2(^+5) 



and is negative in the interval {k'^ , k? 
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